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Abstract: M = 2 supergravity in four dimensions, or equivalently M = 1 supergravity in 
five dimensions, lias an interesting set of BPS solutions that each correspond to a number 
ri^ of charged centers. This set contains black holes, black rings and their bound states, as 

well as many smooth solutions. Moduli spaces of such solutions carry a natural symplectic 
form which we determine, and which allows us to study their quantization. By counting 
\^ the resulting wavefunctions we come to an independent derivation of some of the wall- 

rC^ crossing formulae. Knowledge of the explicit form of these wavefunctions allows us to find 

Tj" quantum resolutions to some apparent classical paradoxes such as solutions with barely 

f — . bound centers and those with an infinitely deep throat. We show that quantum effects 

^^ seem to cap off the throat at a finite depth and we give an estimate for the corresponding 

f^ mass gap in the dual CFT. This is an interesting example of a system where quantum 

^ effects cannot be neglected at macroscopic scales even though the curvature is everywhere 

'k> small. 
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1. Introduction 

One way to try to understand black holes in string theory is to follow them as we slowly 
decrease the value of the string coupling constant. Eventually gravitational interactions 
will become so weak that the black hole will turn into weakly interacting strings and 
branes. The number of states in the latter system with given quantum numbers can often 
be determined via a direct computation in weakly coupled string or field theory. In suitable 
supersymmetric situations, where the number of states can be argued to be independent 
of Qs , this provides a microscopic computation of the entropy of the black hole as was first 
accomplished in [1]. 

In some cases, notably for half-BPS states in AA = 4 SYM [2] [3] and for half-BPS 
states in the D1-D5 system [4] one can do much better and show that the phase space 
of classical smooth supergravity solutions with the same quantum numbers as the black 
hole agrees with (the supersymmetric sector of) the phase space of the weakly coupled 
dual field theory. In particular, the symplectic form obtained in supergravity by restricting 
the full symplectic form to the relevant family of solutions is identical to that obtained 
in the dual field theory [5] [6] [7] [8]. In these cases the entire phase space is subject to 
a non-renormalization theorem and one can see the individual microstates arising from 
supergravity directly. Furthermore, the black hole is to be viewed as a coarse grained, 
thermodynamical description of the corresponding family of smooth solutions [9]. This 
latter statement is the essence of the fuzzball proposal [10, 11, 12, 13]. With such a 
detailed matching between supergravity solutions and the dual field theory at hand we 
can, in principle, ask very detailed questions about black hole geometries. In particular, 
we can try to analyze how and where (semi)classical supergravity breaks down. A precise 
answer to the latter question is bound to shed considerable light on the nature of quantum 
gravity and the information loss paradox. 

Unfortunately, so far full agreement has only been achieved in cases where no large 
macroscopic black holes exist. Motivated by this we consider, in this paper, families of 
solutions of the supergravity equations of motion in four and five dimensions with fixed 
asymptotics that preserve four out of eight supercharges. Our analysis applies both to the 
case of asymptotically flat and asymptotically AdSs solutions. Among such solutions there 
exist large, supersymmetric and extremal black holes, bound states of black holes and black 
rings, as well as many smooth solutions. For a fixed set of given charged centers there is 
a multi-parameter space of BPS solutions that we refer to as the solution space of that 
configuration. At small values of gs such a configuration of charged centers can be related 
to a 0+1 dimensional quiver gauge theory [14], whose BPS moduli space coincides with 
the corresponding solution space in supergravity (i.e. they are parameterized by the same 
variables satisfying the same constraints). This gauge theory describes branes wrapped on 



internal cycles of the compactification manifold and a non-renormalization theorem relates 
the symplectic form derived in this context to the symplectic form on the supergravity 
solution space (i.e. which makes the latter a phase space). Thus, rather than directly 
evaluating the symplectic form in supergravity, we will evaluate it in a weakly coupled 
dual description and rely on the non-renormalization theorem to extend our results to 
strong coupling. We should note that extending the symplectic form in this way would 
not make any sense if the moduli space of solutions was itself not protected by the same 
renormalization theorem. The fact that this space can actually be explicitly computed in 
the small and large Qs regime and matches perfectly is an encouraging confirmation of the 
validity of our approach. 

Obviously, the BPS phase space so obtained is a restricted subset of the full phase 
space of the supergravity theory so great care must be taken in using the associated wave- 
functions. One immediate subtlety that emerges, for instance, is that quantization does not 
completely commute with reduction. Specifically, first imposing a BPS condition and then 
quantizing should give rise to slight deformations of states obtained by first quantizing and 
then imposing a BPS condition. In the latter case the wave functions will also have some 
small support away from supersymmetric solutions [14], whereas in the former they do not. 
Thus the detailed form of the wave-functions will be different. The wave- functions with 
support on the BPS phase space will not be appropriate for many dynamical questions, but 
only for those for which the path- integral localizes on supersymmetric field configurations. 
A more detailed discussion of the subtleties emerging in such a quantization can be found 
in [6]. 

These subtleties notwithstanding, we will find that explicit knowledge of the phase 
space and corresponding wave functions in supergravity is very useful in understanding and 
clarifying the nature of the various supergravity solutions. In particular, for a specific set of 
configurations we will be able to give an alternative derivation of the wall-crossing formula 
which expresses how many states appear or disappear from the spectrum whenever we 
cross a wall of marginal stability [15], including non-primitive crossings. What is more, our 
derivation relies on the explicit phase space structure rather than the ability to decompose 
configurations by tuning moduli at infinity (as in [15]) so we are able to compute the 
number of states for scaling solutions [16] [15] [17] which exist for all values of the moduli 
and hence have no known split attractor flow description. Extending our results to scaling 
solutions with an arbitrary number of centers may shed some light on the idea that these 
solutions are candidate microstate geometries for macroscopic J\f = 2 black holes. This is 
something we hope to come back to in the near future. 

Our results also point to what might well be more general lessons regarding the quan- 
tum structure of black holes. Of particular importance is understanding when quantum 
effects can no longer be neglected. A general lesson that could be drawn from our analysis 
is that the combination of gravity and quantum mechanics can lead to large scale quantum 
structure [18] [19]. We find, for instance, that a small phase space volume (one plank unit) 
can sometimes be spread over a macroscopically large physical volume. There exist a class 
of nominally "classical" solutions with large actions that differ from each other on macro- 
scopic scales but that nonetheless occupy the same plank volume of phase space and thus 



do not represent "good" classical geometries because Heisenberg's uncertainty principle im- 
plies that it is impossible to localize a wavefunction on one of these solutions. This applies, 
in particular, to the aforementioned scaling solutions and implies that, although classically 
they seem to correspond to smooth geometries with infinitely deep throats, quantum effects 
actually cap the throats off at some finite distance. 

By analysing the explicit wavefunctions we discover that the probability density van- 
ishes at the point in phase space that corresponds, classically, to the infinitely deep throat. 
So it seems that, taking into account quantum effects, the infinitely deep throat solution 
is not a good classical solution and, at some depth, the naive classical geometry breaks 
down. We study this using an "effective" capped-throat whose depth is computed in the 
state that localizes as close to the scaling point as possible. In this state we compute the 
expectation value of the size of the "cap" in terms of the charges and relate it to the mass 
gap of the dual CFT. Indeed, without this quantum effect that cuts off the throat, we 
would be forced to understand the strange appearance of a continuous spectrum in the 
CFT because the energy of excitations localized deep down the throat are red-shifted by 
an arbitrary amount when measured from infinity. 

One point to stress is that we will be completely explicit in describing the phase 
space and the corresponding wave functions in cases with two or three centers and also in a 
restricted class of solutions with an arbitrarily large number of centers. The main technical 
tool that allows us full analytic control over these wavefunctions is the observation that, in 
many cases, the solution spaces of multicentered solutions are toric manifolds. Using some 
technology from toric Kahler geometry allows us to perform the geometric quantization of 
these solution spaces explicitly. 

Using some results of [20] our analysis carries over straightforwardly to solutions that 
are asymptotically AdSsxS^xCY. Quantization of such families of solutions yields states 
dual to those in A'^ = (0, 4) MSW conformal field theory, and we hope that this will help 
shed some light on the mysterious nature of this SCFT [21] [22] [20]. 

The outline of this paper is as follows. In section 2, we review the multicenter solutions 
in asymptotically flat and AdSs space-times that we use in the remainder of the paper and 
we summarize the split attractor flow conjecture which allows us to distinguish bona fide 
solutions from ill-defined ones. In section 3 we discuss various aspects of solutions obtained 
in this way with particular emphasis on two-centered and three-centered ones. Section 4 
contains the main result of the paper, which is an explicit expression for the symplectic form 
which is then used to quantize the space of two- and three-centered solutions using methods 
of toric Kahler geometry. We also discuss how to properly incorporate the fermionic degrees 
of freedom in the phase space description and show that these are crucial in order to obtain 
agreement with the wall-crossing formula. In sections 5 and 6 we generalize our procedure 
to certain classes of solutions with an arbitrary number of centers but only up to three 
different types of charges. We again compare to wall-crossing formulae, this time for non- 
primitive charges, and again find perfect agreement. In section 7 we consider what happens 
when centers can move off to infinity in order to clarify which solutions correspond to bound 
and which to unbound states. Finally, in section 8 we discuss the quantum structure of 
solution spaces with a scaling point and compute the way in which quantum effects cap 



off tlie tliroat of these solutions that classically can become infinitely deep. We use these 
results to give some estimate of the mass gap in the dual CFT. We conclude in section 9 
with some comments about the status and meaning of the fuzzball proposal in this context 
and some indications for future work. 

There are also three appendices where we have collected some technical results and 
background material. In appendix A we describe three center solution spaces in full gener- 
ality and detail, proving that the size of the angular momentum is a smooth coordinate on 
that space. In the second appendix we review some results in symplectic toric geometry 
that will be of use in the main text. Appendix C contains a toy model computation of the 
mass gap using AdS/CFT by considering a free scalar field on an approximate geometry. 

2. Review of Multicentered Black Holes 

Let us begin with a brief review of multicentered black hole solutions oiN = 2 supergravity 
in 4 dimensions. From the 10 dimensional point of view we will be considering Type IIA 
on a Calabi-Yau threefold X (though most of our results carry over to IIB). Furthermore 
we shall work in the limit where the volume of X is very large so we can restrict ourselves 
to just the cubic part of the prepotential. 

The multicenter solutions are determined by specifying a number of charges, Fq, and 
distributing their location in the spatial M^. These charged centers correspond in the 
10 dimensional picture to branes wrapping even cycles in the CY3. There are 262 + 2 
independent such cycles in homology, with 62 the second Betti- number of X, each giving 
rise to one of the 262+2 abelian vector fields of the N = 2 supergravity. We will often denote 
brane configurations by their coefficients in a basis of homology, i.e. F = {p^ ,p , qa: Qo) = 
pO -\.pA£)^ _|_ q^f)^ -\- qqV, where the Da form a basis of ii^{X, Z), the Da make up a dual 
basis and V is the unit volume element of X, dual to 1. There is a natural intersection 
product on H*(X, Z) given by 

(Fi,F2) := -p'iqi+ptqA " qAP2 + qIpI (2-1) 

Furthermore the moduli of the Calabi-Yau appear as scalar fields in the four-dimensional 
effective theory. In the solutions we will be considering the hyper-multiplet moduli will be 
kept constant while the ones in the vector-multiplet will be allowed to vary dynamically. 
It is thus necessary to choose boundary values for the scalar moduli at infinity. We will 
briefly discuss the split attractor flow conjecture which relates the existence of solutions at 
particular values of the moduli at infinity to the existence of certain flow trees in moduli 
space. 

We will be very brief in our review as these solutions are discussed in great detail 
in the references [23], [24] and [15]. Although below we will discuss primarily the four 
dimensional version of these solutions it is well known that these solutions can be uplifted 
to five dimensions [25] yielding multicentered solutions with a U{1) isometry and a Gibbons- 
Hawking base that can be made smooth by required all the four-dimensional charges to 
come from fluxed Z)6-branes [26] [27]. These are solutions of 11 dimensional supergravity 
compactified on the same Calabi-Yau, X. 



2.1 Four Dimensional Solutions 

Our starting point is the class of multicentered black holes solutions first discussed in 
[28] [23] and later explicitly constructed in [24]. They are given by the following four 
dimensional metric, gauge fields and moduli: 

ds^ = -^{dt + ojf + j:dx'dx\ 

A'> = ^^{dt + uj)+uo, (2.2) 

.A (9 log S , , , , 
A"^ = -—^ {dt + uj)+Ad, 
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All functions appearing in this solution can be expressed in terms of 262 + 2 harmonic 
functions: 

H = {H'>,H^,Ha,Ho) = V , ^ - 2Im(e-^°[7)U , (2.3) 

a ' ' 

where the components of the Ta take values in H*^(X,'Ij), the integral even cohomology of 
the Calabi-Yau X, e*" is the phase of the total central charge at infinity (i.e. ^(r)|c 

_Z_\ 

\z\\ 



{Y^s^sMoo) and e*- = ^|oo) and ^ = -^^. 



The function S appearing in the metric in (2.2) is known as the entropy function. 
When evaluated at Xa it is proportional to the entropy of the black hole whose horizon lies 
at Xa- This follows from the Bekenstein-Hawking relation and the fact that S determines 
the area of the horizon at Xa- 

The rest of the solution is defined in terms of the harmonic functions as follows 

dujQ = -kdH^ , 
dAJl = *dH^ , 
*dLO = {dH, H) 

If we take the prepotential to be cubic, which is tantamount to taking the large volume 
limit in IIA, we can find an explicit form for the entropy function [29] 



(2.4) 



V (^0)2 ' 

Q^ = (^I?ABcy^V)^ 

DABcy^y'' = -2HcH'^ + DabcH^H'' , 

where one should note that the Hodge star * is the one of flat M'^. Finally there are A^ — 1 
consistency conditions on the relative positions of the A^ centers. These conditions come 



from integr ability of the equations of motion for cu. They take the simple form 

(^,r,)U=,, = 0, (2.5) 

or written out more explicitly^ 

Y, ^^^^^^ = {h,T,) , (2.6) 

a, a^b 

where rab = \xab\ = \xa - Xb\. 

An important property of a configuration with a sufficient number of centers is that 
although the centers bind to each other there is some freedom left to change their relative 
positions. These possible movements can be thought of as flat directions in the interaction 
potential. Equation (2.6) constrains the locations of the centers to the points where this 
potential is zero. As, for a system with A'^ centers, there are A — 1 such equations for 3 A — 3 
coordinate variables (neglecting the overall center of mass coordinate) there is, in general, 
a 2A — 2 dimensional moduli space of solutions for fixed charges and asymptotics. This 
space may or may not be connected and it may even have interesting topology. We will 
refer to this as the moduli space of solutions or solution space; the latter terminology will 
be preferred as it is less likely to be confused with the moduli space of the Calabi-Yau in 
which the scalar fields t take value. The shape of this solution space does, in fact, depend 
quite sensitively on where the moduli at infinity, i |ooi he in the Calabi-Yau moduli space 
(as the latter determine h on the RHS of eqn. (2.6)). Quantization of some low-dimensional 
instances of this space will be the primary goal of the rest of the paper. 

Let us briefiy recall, once more, that these four-dimensional solutions can be up- 
lifted to five-dimensions [25]. Moreover, for certain choices of charges, namely pure D6's 
with only abelian flux turned on^, the five-dimensional uplift of the solutions is smooth 
[30] [27] [31] [32]; there are, at worst, conical singularities at the location of each center. The 
five dimensional solutions are, in general, topologically non-trivial multi-Taub-NUT spaces 
with fiuxes supported on the two cycles between the Taub-NUT centers. The fact that such 
smooth solutions, with neither horizons nor singularities, can be found is important for the 
fuzzball proposal and some of the discussion that follows will involve such configurations 
within the context of this proposal. For a more details and discussion of these smooth 
solutions the reader is referred to the references given above. 

2.2 Angular Momentum 

These solutions are stationary and, as shown in [23], they carry an angular momentum 

^In this sum we imply summation over a different from some fixed b. In case we don't give a separate 
mention of tlie summation index we mean a double sum. E.g. '}2a<b ^^ ^^'^ over a and b with a smaller 
than b, while "^^ ^^^ would mean a sum over only a such that a is smaller than some fixed b. 

^Note that in the duality frame we are working in it is purely fluxed D6-centers that give rise to smooth 
five dimensional solutions. In other duality frames this can be purely fluxed D4's, etc. It is probably more 
precise to consider zero entropy centers instead of only smooth centers, as the former is a duality invariant 
statement. 



equal to 

j^lj2{Ta,n)Xab_ (2.7) 

^ a<b """-^ 

Note that J is quantized in half integer units. 

As angular momentum will play an important role in the rest of this paper, providing a 
natural coordinate on the solution space, we will derive some more useful ways of expressing 
it. Multiplying the condition (2.6) by Xa and then summing over the different centers shows 
that 

J=\Y.'^KVa)xa. (2.8) 

a 

By using the fact that X]a(/i,ra) = one can rewrite the last expression as 

J=\ 5^(/i,r,)x,6. (2.9) 

Starting from this formula we can show that the size of the angular momentum can be 
compactly written in terms of the inter-center distances Tab- Squaring (2.9) gives 

J'' = \Y. ^ ^^' ^») ^^' r-) ^-^ ■ ""-^ ■ (2.10) 



Now we can use that for any three points labeled by a, b, c there is the relation Xab ■ x^b 

\{^lb + ^lb-^lc) 



^('"afe + '"cfe ~ ''^'ac) ^^d a little more algebra reveals that 



\J\ = 2J-^^^'^-)^^^'^b)^lb- (2-11) 

y a<b 

2.3 The Split Attractor Flow Conjecture and Wall Crossing 

So far we have reviewed a class of four dimensional solutions but these solutions are rela- 
tively complicated and it is non-trivial to determine if they are well-defined everywhere. In 
particular the entropy function, S, that appears in the solution involves a square root and 
may take imaginary values in some regions (when uplifted to five dimensions this can lead 
to closed timelike curves [31] [27] [26]). In [23] and [15] a simplified criterion was proposed 
for the well behavedness of such solutions which we will now relate (very briefly) . 

In [23] a conjecture was proposed whereby the existence of multicentered solutions is 
equivalent to the existence of an attractor flow tree. The latter is a graph in the moduli 
space beginning at the moduli at infinity and ending at the attractor points for each 
center. The edges correspond to single center flows towards the attractor point for the sum 
of charges further down the tree. Vertices can only occur where single center flows (for a 
charge F = Ti -|- r2) cross walls of marginal stabilities so that the central charges are all 
aligned {\Z{T)\ = \Z{Ti)\ + \Z{T2)\). The actual (spatially dependent) flow of the moduU, 
t , for a multicentered solution will then be a thickening of this graph. A sample tree is 
displayed in figure 1. 

As mentioned, the main purpose of the attractor flow tree is to allow us to determine 
if a solution is well defined. For a single centered black hole the entropy function, E, 



undergoes a monotonic flow from inflnity to the horizon. At infinity the value of S depends 
on the choice of modufi (boundary conditions) while at the horizon it flows to a flxed 
value depending only on the charges, as the moduli are flxed by an attractor mechanism. 
Spherical symmetry dictates that the moduli depend only on a radial variable so the flow 
through moduli space is indeed just a single line from the moduli at inflnity to the attractor 
value. If S should become imaginary somewhere along this flow the solution would suffer 
from pathologies. However, since the flow is monotonic it need only be checked at its initial 
(the moduli at inflnity) and flnal points (the attractor point). 

For a multicentered system the moduli depend on three variables and the flow is no 
longer monotonic in a straightforward way (it is not even a one dimensional tree but rather 
a "fat graph"). By assuming that solutions could be built constructively by bringing in 
centers from inflnity [23] was able to conjecture that even for multicentered conflgurations 
we can study a flow tree in the moduli space (recall the actual flow will be a "fat" version 
of this) and study each leg of the flow to check for pathologies. The conjecture is then 
that if the tree exists (each leg is pathology free) then the full solution is actually well 
behaved (see [23] [33] for more details). There is considerable evidence for this conjecture 
[23] [14] [33] [15] and our computation in Section 4.3.4 will provide even further support. 

Although our treatment of attractor trees here is very brief we will discuss one more 
important result from [15], the so called wall crossing formulae. These are formulas for the 
index of BPS states of given total charge. They are derived using that the index of the 
BPS Hilbert space corresponding to an attractor flow tree can be determined in terms of 
contributions from each vertex and node. A given vertex for a split T -^ Ta + Tb contributes 
a factor |(ra,rfe)|, associated to the number of states in the angular momentum multiplet 
coming from the crossed electric and magnetic flelds of the centers Pq and Pf,. Each node, 
meanwhile, contributes a factor ri(Pa) where $^(Pa) represents the index associated to 
the charge P^ itself (i.e. coming from "internal" degrees of freedom and determined by 
the entropy formula (2.4) with H being replaced by Pq). Thus, for instance, the index 
corresponding to flgure 1 is 

|(P3,Pl + P2)| |(Pl,P2)| l^(Pl) n{T2) 0(P3) (2.12) 

Note that the spacetime contribution (from the vertices) and the "internal" contributions 
(nodes) are easily and clearly separated in this computation. For more details, including a 
derivation (assuming the split attractor conjecture) the reader is referred to [23] and [15]. 

3. Simple Solution Spaces 

Let us describe some simple moduli spaces of solutions in order to have some feeling for 
the spaces we wish to quantize (in section 4). We begin with the simple case of the two 
centers solution and then discuss the three centers case. 

3.1 The Two Center Case 

The solution space for two centers, when it exists, is two dimensional. The constraint 




Figure 1: Three centered attractor flow tree. The system is composed of three centers of charge 
Fi, r2 and T^ and the moduH at inflnity are at the value labelled by the yellow circle. Each leg of 
the tree above represents a single center flow towards the attractor value associated with the total 
charge below that point. Thus the first flow is towards the attractor point for F = Fi + r2 + Fa. 
After the first split the flows are towards the attractor points for charges F3 and F4 = Fi + F2. In 
each case the split occurs along walls of marginal stability (thick blue lines). The flrst, horizon, line 
of MS corresponds to \Z{T)\ = |Z(F3)| + |Z(F4)| while the second is for \Z{Ti)\ = |Z(Fi)| + |Z(F2)|. 



equation (2.6), 

(ri,r2 



2;i2 



(/i,ri) (3.1) 



admits solutions if {Ti,T2)/{h,Ti) is positive and finite. These solutions have fixed inter- 
center distance, X12, and are hence parameterized only by the direction of the inter-center 
axis (recall that we are factoring out the center-of-mass coordinate). Thus the solution 
space is topologically an 5^. 

If (ri,r2) = {h,Ti) = then the centers are unbound and can be located anywhere. 
In this case the solutions carry no angular momentum and are static. To quantize static 
solutions one needs to include velocities, as was e.g. done in [34]. Whether or not static 
solutions give rise to proper BPS states is an interesting and subtle question on which we 
further comment in our conclusions in section 9. 



3.2 The Three Center Case 

The three center solution space is four dimensional. Placing one center at the origin (fixing 
the translational degrees of freedom) leaves six coordinate degrees of freedom but these 
are constrained by two equations. This leaves four degrees of freedom, of which three 
correspond to rotations in S0(3) and one of which is related to the separation of the 
centers. 
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= -{hJ:, 



The constraint equations take the form 

=: a (3.2) 

=: (3 (3.3) 

in a self-evident notation. The nature of the solution space simplifies considerably if either 
aoT f3 vanish so let us first consider this case (if both vanish there is an overall scaling degree 
of freedom and the centers are unbound; we will, in section 7, see that this corresponds to 
a degenerate symplectic form and is thus not amenable to quantization using the methods 
described in this paper). For definiteness we will take a = 0; in this case Ylpi^^^p) — 
which implies {h,T2) = [i. Thus from (2.8) we find 

/=^r23S (3.4) 

with z defined to be parallel to X2 — X3. 

The solution has an angular momentum vector J* directed between the centers 2 and 3 
and the direction of this vector defines an S^ in the phase space which we will coordinatize 
using 9 and (p. The third center is free to rotate around the axis defined by this vector (since 
this does not change any of the inter-center distances) providing an additional C/(l), which 
we will coordinatize by an angle o", fibred non-trivially over the S^. Finally the angular 
momentum has a length which may be bounded from both below and above and this 
provides the final coordinate in the phase space, j = \J\. This construction is perhaps not 
the most obvious one from a spacetime perspective but, as we will see, in these coordinates 
the symplectic form takes a simple and convenient form. When a = it is clear from (3.4) 
that j is a good coordinate on the solution space but this is not immediately obvious for 
the more complicated case of a 7^ 0. This is nonetheless true and, as shown in Appendix 
A, this is always a good coordinatization of the three center solution space (though for 
a 7^ the relation between (j, a, 9, (p) and the coordinates Xp is not as straightforward). 

The quantization of these solutions is particularly interesting in certain cases where 
classical reasoning leads to pathologies. Before proceeding with quantization we will first 
briefly review these solutions. 

3.2.1 Scaling Solutions 

As noted in [35] and [15], for certain choices of charges it is possible to have points in the 
solution space where the coordinate distances between the centers goes to zero. Moreover, 
this occurs for any choice of moduli so it is, in fact, a property of the charges alone. 

Such solutions occur as follows. We take the inter-center distances to be given by 
i"ab = AFaf, -|- 0{X'^) (fixing the order of the ab indices by requiring the leading term to 
always be positive). As A —> we can always solve (3.2) by tweaking the A^ and higher 
terms. The leading behaviour will be r^f, ~ XTab but clearly this is only possible if the 
Tab satisfy the triangle inequality. Thus any three centers with intersection products Tab 
satisfying the triangle inequalities define a scaling solution. 
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We will in general refer to such solutions as scaling solutions meaning, in particular, 
supergravity solutions corresponding to A ~ 0. The space of supergravity solutions con- 
tinuously connected (by varying the Xp continuously) to such solutions will be referred to 
as scaling solution spaces. We will, however, occasionally lapse and use the term scaling 
solution to refer to the entire solution space connected to a scaling solution. We hope the 
reader will be able to determine, from the context, whether a specific supergravity solution 
or an entire solution space is intended. 

These scaling solutions are interesting because (a) they exist for all values of the 
moduli; (b) the coordinate distances between the centers go to zero; and (c) an infinite 
throat forms as the scale factor in the metric blows up as A^^. Combining (b) and (c) we see 
that, although the centers naively collapse on top of each other, the actual metric distance 
between them remains finite in the A — > limit. In this limit an infinite throat develops 
looking much like the throat of a single center black hole with the same charge as the total 
charge of all the centers. Moreover, as this configuration exists at any value of the moduli, 
it looks a lot more like a single center black hole (when the latter exists) than generic non- 
scaling solutions. As a consequence of the moduli independence of these solutions it is not 
clear how to understand them in the context of attractor flows; the techniques we develop 
in this paper provide an alternative method to quantize these solutions that applies even 
when the attractor tree does not. 

Unlike the throat of a normal single center black hole the bottom of the scaling throat 
has non-trivial structure. If the charges, Fq, are zero entropy bits (e.g. D6's with abelian 
flux) then the 5-dimensional uplifts of these solutions will yield smooth solutions in some 
duality frame and the throat will not end in a horizon but will be everywhere smooth, even 
at the bottom of the throat. Outside the throat, however, such solutions are essentially 
indistinguishable from single center black holes. Thus such solutions have been argued to 
be ideal candidate "microstate geometries"^ corresponding to single center black holes. In 
[15] it was noticed that some of these configurations, when studied in the Higgs branch of 
the associated quiver gauge theory, enjoy an exponential growth in the number of states 
unlike their non-scaling cousins which have only polynomial growth in the charges. 

3.2.2 Barely Bound Centers 

We will be brief here as such configurations are discussed in more detail in section 7. For 
certain values of charges and moduli we will see that it is possible for some centers to 
move off to infinity. Although this would normal signal the decay of any associated states 
(as happens, for instance, for two centers at a wall of marginal stability [23]) we will see 
that this is not always the case. In particular, it is important to distinguish between cases 
when centers are forced to infinity (marginal stability) versus those where there is simply 
an infinite (flat) direction in the solution space (threshold stability; see [20, Appendix 
B]). Although the first case clearly signals the decay of a state, in the second case, when 
centers move off to infinity along one direction of the solution space but may also stay 



^The unfortunate terminology "microstate geometry" is not intended to suggest that the classical solution 
necessarily corresponds to a black hole microstate. Rather, as we will see later, quantizing the space of such 
solutions as a phase space provides a set of quantum states that may or may not have a classical realization. 
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within a finite distance in other regions of the solution space, it is still possible to have 
bound states, as we will demonstrate. Quite essential to this argument is the fact that in 
some cases, although the solution space may seem naively non-compact (in the standard 
metric on M?^^"^), its symplectic volume is actually finite and it admits normalizable wave- 
functions (whose expectation values, we argue in section 7, are finite). There are also cases 
with unbound centers where the symplectic form on the solution space is degenerate and, 
in such cases, it is not clear if there is a bound state. Such cases are not amenable to 
treatment by the methods developed here, see also the discussion on this point in section 



4. Symplectic Structure and Solution Space Quantization. 

In this section we study the moduli space of our solutions using a symplectic form which 
is derived from the quiver quantum mechanics action of multiple intersecting branes in the 
coulomb branch [14]. We determine this symplectic form explicitly for the three center case 
and show it to be non-degenerate (with some important exceptions where degenerations 
have a clear physical interpretation) and proceed to count the number of states in this 
moduli space via geometric quantization. The estimated number of states at large charge 
nicely matches the expected number of states from the wall-crossing formula given in [15]. 

Motivated by a non-renormalization theorem [14] and the exact agreement of our state 
counting with Denef and Moore's wall-crossing formula we propose that the same symplec- 
tic from should be derivable from the super-gravity action following the logic in [5] (see also 
[36] and references therein). Actually a similar structure to the quiver quantum mechanics 
symplectic form emerges from the gauge field contribution of the supergravity action. So 
we can rephrase our conjecture in the following way: the other putative terms contributing 
to the symplectic form from supergravity cancel or only change the normalization '^. 

Note that it is not a priori clear that our moduli space of solutions corresponds to 
a full phase space rather than a configuration space. The fact that the solution spaces 
are always even, (2N — 2), dimensional is a positive indication that they are indeed phase 
spaces as is the fact that the symplectic form is non-degenerate in the two and three center 
cases. Moreover, when degenerations do occur they have well-defined physical meaning 
further supporting this claim. 

4.1 Symplectic Form from Quiver QM 

In order to be able to quantize the space of solutions described above, we need to determine 
its symplectic form which is induced from the supergravity Lagrangian. The idea is very 
simple, one first determines the symplectic structure on the full space of solutions to the 
supergravity equations of motion, which takes the form 

^= /^S/'^f^J^V'^'A^- (4-1) 



d{di^^) 



*The last possibility is seen for example in [5]. 
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where L is the relevant supergravity Lagrangian. Here, the integral is over a Cauchy surface, 
(f) represents a basis of the fields that appear in the Lagrangian and it is assumed that 
the Lagrangian does not contain second and higher order time derivatives. The induced 
symplectic form on the space of BPS solutions is then simply the restriction of (4.1) to 
that space. 

In general, there is no guarantee that the symplectic form obtained in this way is non- 
degenerate. If it turns out to be degenerate then we should have included further degrees 
of freedom to arrive at a non-degenerate symplectic form. This happens for example when 
all centers are mutually BPS, i.e. if all inner products (Tp^Vq) vanishes. In this case the 
symplectic form, as we will see below, is identically zero, and the solution space is therefore 
better thought of as being a configuration space. In order to obtain a non-degenerate 
symplectic form, we could try to include, for example, small velocities for the centers [34]. 
It is not clear, however, whether this can be done while saturating the BPS bound. We 
return to this point briefly in our concluding remarks in section 9. 

It may sound surprising that generically the BPS solutions spaces will carry a non- 
degenerate symplectic form, since they are all time-independent. Crucially, the solutions we 
consider are stationary but in general not static. Stationary solutions do carry non-trivial 
momentum despite being time-independent and this is what gives rise to the non-trivial 
symplectic form. As a consequence of this we will see that the Hilbert space decomposes 
into angular momentum multiplets. The symplectic form (4.1), when evaluated on a family 
of static solutions, will simply vanish. 

The idea to quantize spaces of BPS supergravity solutions using the restriction of (4.1) 
was successfully applied in [5] [6] [8] (see also [37] for an extensive list of older references and 
[38] [39] [36] for more recent, similar work). However, when we try to apply similar methods 
to the space of multicentered black hole solutions, the expressions become very lengthy 
and tedious and we run into serious technical difficulties due to the complicated nature 
of the multicentered solutions. Therefore we will proceed differently and try to derive the 
required symplectic form from a dual open string description of the supergravity solutions. 

From the open string point of view, which is appropriate for small values of the string 
coupling constant, black hole bound states correspond to supersymmetric vacua of a suit- 
able quiver gauge theory. The connection between the supergravity solutions and gauge 
theory vacua becomes clear once we study the Coulomb branch of the gauge theory. For 
simplicity we will assume that all centers have primitive charges; the extension to non- 
primitive charges is straightforward, as the supergravity solutions are not specifically sen- 
sitive to whether charges are primitive or not. 

For primitive charges, the quiver gauge theory has as gauge group a product of J7(l)'s, 
one for each center, plus bifundamentals. The number of bifundamentals minus the number 
of anti-bifundamentals charged under a pair of C/(l)'s is equal to (Fp, Fg), the inner product 
of the charges. To study the Coulomb branch of this theory, we integrate out massive 
bifundamentals, resulting in an effective theory for the vector multiplets only. The degrees 
of freedom of the vector multiplet in -|- 1 dimension consist of three scalars x*, four 
fermions, one auxiliary field D and one gauge field A. As was shown in [14], the low-energy 
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effective action for A^ abelian vector multiplets necessarily takes the form 

L = y {—UpDp + Ap ■ Xp) + fermions + higher order terms, (4-2) 

p 

where p labels a brane localized at Xp and 

VpC/g = VgUp = -(Vp xAg + VgX Ap). (4.3) 

This form of the action is determined by requiring that it is invariant under supersymmetry. 
The solutions to (4.3) are given in terms of harmonic functions, in particular we find that 

Up = {Tp,H{xp)) (4.4) 

where H = + ^ , _^ , and the pole at x = Xp in H does not contribute to (4.4). 
We therefore see that except for the tree-level constant term 0, there are only one-loop 
contributions to Up and Ap. Supersymmetry prohibits higher loop contributions to Up and 
Ap. Luckily, we only need Up and Ap to match the gauge theory to supergravity and to 
extract the symplectic form, and since they do not receive higher loop corrections we can 
safely extrapolate the results from the gauge theory regime at small gs to the supergravity 
regime at large Qg. 

It would be interesting to rephrase the above non-renormalization theorem more di- 
rectly in superspace. The superspace formulation of the -|- 1-dimensional vector multiplet 
is in terms of a linear superfield S^^ with D^S = Z)^S = [40] [41] . It is therefore tempting 
to try to write (4.2) in terms of a twisted superpotential 

L= I dO'^de^Hc.pi^) (4.5) 

with H some sort of prepotential for Up and Ap. A first investigation indicates that this is 
not possible; if it were possible, then it would be clear that the string coupling, which can 
be viewed as the expectation value of the lowest component of an auxiliary spurious chiral 
superfield, cannot contribute to the twisted superpotential. Nevertheless, one can verify 
explicitly that it is not possible to write down first order terms like (4.2) which involve a 
non-trivial coupling between a linear superfield and a chiral superfield, thereby giving an 
alternative proof of the non-renormalization theorem. 

Given (4.2), the space of supersymmetric vacua of the gauge theory is given by the so- 
lutions of the D-term equations f/p = 0. These are identical to the supergravity constraint 
equations (2.5), (2.6), as 9p can be identified with the constant term in the supergravity har- 
monic functions. The symplectic form then follows immediately by restricting ^ 5xp/\5Ap 
obtained from (4.2) to the solution space C\p{Up = 0}. This symplectic form is independent 
of Qs and must therefore agree with the symplectic form obtained from supergravity. 

Phrased in terms of the supergravity solution, the symplectic form becomes 

^ = \Y.K^^^P^^^d{xp))- (4.6) 
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where Ad is the "spatial" part of the 4d gauge field defined in (2.2). Notice that 

SAliixp) = (5^j;(x))U=,^ + {6x'pdkAi,{x))\^=.^. (4.7) 

To proceed further we denote 

fp = r^-^- (4.8) 

"V" — "V" 

With some work, one can show that 

6Ad = Y, TpeijkSx; djfp dx'' (4.9) 

p 

Using the above form of 6 Ad and the fact that d'^A^^ix) can be replaced by ^{d^Ali{x) — 
d^A^{x)) = ^J^ll^{x) we finally obtain for the symplectic form 

0_1V^/T T^ \ (^ijk{S{Xp - XgY A 5{Xp - XgY) {Xp - Xqf 

^ pH '^ ^' 

Note that the overall translational mode does not contribute to this symplectic form. 

The symplectic form (4.10) applies for any number of centers but it must still be 
restricted to the solution space defined by the constraint eqns. (2.6). Since these spaces 
are quite complicated we will only be able to analyse this restriction for the two and three 
center case and some simple examples with arbitrary number of centers. 

4.2 The Solution Space as a Phase Space 

The symplectic form of the phase space must be closed and non-degenerate in order for 
the latter to be a well-defined symplectic manifold. We would like to demonstrate these 
properties for Cl in general but the second property depends on the structure of the solution 
space (and in fact does not hold in some degenerate cases as we will see) so we will only 
be able to demonstrate it for the two and three center cases. As previously mentioned it is 
not a priori evident that the solution space is a phase space so in principle Jl might have 
been degenerate. 

A direct calculation shows that dCt is closed as a two-form on the solution space 

dn^Y. ^T^ f ^''' - 3^^f^^) 8x],g A 6x1^ A 5xlg = (4.11) 

p^g \^pq\ \ Fpgl / 

The last equality is most easily worked out in a coordinate basis. 

4.3 Quantization of the Solution Space 

In sections 3.1 and Appendix A we described the solution space of two and three-centered 
solutions. Here, we will described the quantization of these solution spaces using the 
symplectic form (4.10). It turns out that the structure of the solution space is determined 
to a large extent by the symmetries of the problem. Recall that we already removed the 
overall translational degree of freedom from the solution space, which we can do for example 
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by fixing one center to be at the origin, e.g. xi = 0, or by fixing the center of mass to 
be at the origin. The translational degrees of freedom and their dual momenta give give 
rise to a continuum of BPS states, but this continuum yields a fixed overall contribution 
to the BPS partition function. By factoring out this piece we are left with the reduced 
BPS partition function, and the quantization of the solution space we consider here yields 
contributions to this reduced BPS partition function. 

Besides the overall translational symmetries, the constraint equations (2.6) are also 
invariant under global 50(3) rotations of the Xp. These rotations do appear in a non- 
trivial way in the symplectic form. Indeed, if we insert 6xp = e*°''n"Xp, which corresponds 
to an infinitesimal rotation around the n-axis^, in (l we obtain 

Q -L ^— ^ , ^^ijk^iab"^ \^p ^qj 0\Xp Xq) [Xp Xqj 

~ 4 2^\ P' 9/ 1^ -Y |3 

= v}5f (4.12) 

where J* are the components of the angular momentum vector (see section (2.2)) 

•^^ = lE(rp,r,)^^ (4.13) 

The fact that the symplectic form takes the simple form in (4.12) is not surprising. This 
is merely a reflection of the fact that angular momentum is the generator of rotations. 
Fortunately, this simple form completely determines the symplectic form in the two and 
three-center case. 

4.3.1 Two-Center Case 

The two-center case is easy to describe. There is only a regular boundstate for (ri,r2) 7^ 
and {h,Tp) 7^ 0, and the constraint equations immediately tell us that X12 is fixed and 
given by 

In other words, xi — X2 is a vector of fixed length but its direction is not constrained. The 
solution space is simply the two-sphere, and the symplectic form is proportional to the 
standard volume form on the two-sphere. In terms of standard spherical coordinates it is 
given by 

n= ^{Ti,T2)sm9d9Ad(t)= \J\ sin 9 d9 A dcf). (4.15) 

We can now quantize the solution space using the standard rules of geometric quantization 
(see e.g. [42] [43]). We introduce a complex variable z by 

z^ = '-±^e'^^ (4.16) 

1 — cos 9 



^In other words, we compute tx^ with X the vector field ^ e^°'''n'^Xp-^ 
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and find that the Kahler potential corresponding to 0, is given by 
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/C = -2|J|log(sing) = -|J|log K^^^_^J . (4.17) 

The holomorphic coordinate z represents a section of the line-bundle C (over S^, the so- 
lution space) whose first Chern class equals r2/(27r). The Hilbert space consists of global 
holomorphic sections of this line bundle and a basis of these is given by ipm{z) = z™. How- 
ever, not all of these functions are globally well-behaved. For example, regularity at z = 
requires that m > 0, and to examine regularity at z = oo one could e.g. change coordinates 
z -^ 1/z and use the transition functions of £ to find out that m < 2\J\. Equivalently, we 
can directly examine the norm of -i/'m by computing 

\lPm\^r. f dY0\e-^\^m{z)\'' (4.18) 



where dvol is the volume form induced by the symplectic form. In our case we therefore 
find 

\^m? ^ / rfcos^#(l + cos0)l^l+'"(l-cose)l-^l-™ (4.19) 

and clearly f/^m only has a finite norm if — | J| < m, < | J|. The total number of states equals 
2| J| -|- 1. This is in agreement with the wall-crossing formula up to a shift by 1. We will 
see later that the inclusion of fermions will get rid of this extra constant. 

The integrand in (4.19) is a useful quantity as it is also the phase space density asso- 
ciated to the state V'm- According to the logic in [9, 44] the right bulk description of one 
of the microstates ipm should be given by smearing the gravitational solution against the 
appropriate phase space density, which here is naturally given by the integrand in (4.19). 
We will come back to this point later, but observe, already, that since there are only 2| J| -|- 1 
microstates, we cannot localize the angular momentum arbitrarily sharply on the 5^, but 
it will be spread out over an area of approximately vr/| J| on the unit two-sphere. It is 
therefore only in the limit of large angular momentum that we can trust the description of 
the two-centered solution (with two centers at fixed positions) in supergravity. 

4.3.2 Three Center (Non-Scaling) Case 

We now return to the symplectic form (4.12). In the three-center case, we expect four 
degrees of freedom. As discussed in section 3.2, three of those are related to the possibilities 
to rotate the system, whilst the fourth one can be taken to be the size of the angular 
momentum vector \J\. The specific form of (4.12) strongly suggests that these are also the 
variables in which the symplectic form takes the nicest form. 

We therefore take as our basic variables J* and o", where a represents an angular 
coordinate for rotations around the J-axis. Obviously, a does not correspond to a globally 
well-defined coordinate, but rather should be viewed as a local coordinate on an 5^-bundle 
over the space of allowed angular momenta. Ignoring this fact for now, the rotation 6xp 

4 



^lab^a^b i]^g^^ ^g used in (4.12) corresponds to the vector field 



_ri^d_ijkjjk_^ (A 20) 

"" \J\ da^' -^ dJ^- ^ ^^' 
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The second term is obvious, as J is rotated in the same way as the Xp. The first term 
merely states that there is also a rotation around the J-axis given by the component of n 
in the J-direction. The final result in (4.12) therefore states that 

0(X„,m^--) = nW n(X„,— )=0. (4.21) 

aJ* oa 

It is now easy to determine that 

B B T^ ~ B B P 

^(7^,7^) =ei7fc^, ^{i^,^) = -r^- (4.22) 

Denoting \J\ as j, and parameterize J* in terms of j and standard spherical coordinates 
0, (^, the symplectic form defined by (4.22) becomes 

O = J sin ede Ad(j)- dj A da. (4.23) 

However, we clearly made a mistake since this two-form is not closed. The mistake was 
that a was not a well-defined global coordinate but rather a coordinate on an S^-bundle. 
We can take this into account by including a parallel transport in a when we change J*. 
The result at the end of the day is that the symplectic form is modified to 

n = j sin edeAd(t>- dj A Da (4.24) 

with Da = da — A, and dA = sin 9 d9 A dcj), so that A is a standard monopole one-form on 
5^ . A convenient choice for A is A = — cos 9 dcj) so that finally the symplectic form can be 
written as a manifestly closed two-form as^ 

n = -d{j cos 9) Ad(f)- dj A da. (4.25) 

This answer looks very simple but in order to quantize the solution space, we have to 
understand what the range of the variables is. Since 0, (p are standard spherical coordinates 
on 5^, (j) is a good coordinate but degenerates at 9 = 0,7r. The magnitude of the angular 
momentum vector j is bounded as can be seen from (4.13). By carefully examining the 
various possibilities in the three-center case (see appendix A), one finds that generically j 
takes values in an interval j G [j_ , j+J , where j = j_ or j = j+ only if the three-centers lie 
on a straight line. An exceptional case is if j_ = implying that the three-centers can sit 
arbitrarily close to each other as seen in appendix A. Note that this latter case corresponds 
exactly to the scaling solutions described in section 3.2.1. 

As we mentioned above, at j = j_ and j = j_|_ the centers align, and rotations around 
the J-axis act trivially. In other words, at j = j± the circle parametrized by a degenerates. 
Actually, we have to be quite careful in determining exactly which U{1) degenerates where. 
Fortunately, what we have here is a toric Kahler manifold, with the two U{1) actions given 
by translations in (p and a, and we can use results in theory of toric Kahler manifolds from 
[45] (see also [46] and Appendix B) to describe the quantization of this space. 



^This result is re-derived in a more straightforward but tedious way in Section 6. 
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We start by defining x = j and y = jcos6 to be two coordinates on tlie plane. Then 
the ranges of the variables x and y are given by 



jL > 0, j+-x> 0, 



X 



y > 0, x + y>0. 



(4.26) 



Together these four^ inequalities define a Delzant polytope in M? which completely specifies 
the toric manifold (see Appendix B). At the edges a U{1) degenerates, at the vertices both 
C/(l)'s degenerate. The geometry and quantization of the solution space can be done purely 
in terms of the combinatorial data of the polytope (see figure 2). We have to distinguish 
two cases. First, when j_ = 0, which corresponds to a scaling point inside the solution 
space, our toric manifold is topologically CP?^ ^ g)- I^ the case j_ > 0, where a scaling point 
is absent, the solution space becomes the blow-up of CPq ^ 2)) which can be identified as 
the second Hirzebruch surface F2. These statements can all be verified by e.g. constructing 
the normal fan to the relevant polytopes. 



- j cos 6 




Figure 2: (Left) The polytope for j'_ = corresponding to 



J2 
(1,1.2)- 



(Right) The polytope for 



j_ > 0. This corresponds to the second Hirzebruch surface, F2, a blow-up of ( 



(1,1,2)- 



For the purpose of quantizing the system we will assume that j- > 0; we will return 
to the case of j- = in section 8. Thus the results of the rest of this section only apply for 
j_ > 0. Furthermore we assume that all three centers carry different charges; if two centers 
carry identical charge one needs to take into account their indistinguishability, quantum 
mechanically, and take a quotient of the corresponding solution space. We won't consider 
this possibility in this section but come back to it in detail in section 5. 

The construction of canonical complex coordinates is done by constructing a function 
g as 

1 



9= ^[{x- i_) log(x - J-) + (j+ - x) log{j+ - x) 
+{x - y) log(a; - y) + {x + y) log(x + y)] 



(4.27) 



^In the case j_ — the first equation is redundant and is not part of the characterization of the polytope. 
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which is related in an obvious way to the four inequahties in (4.26). Then the complex 
coordinates can be chosen to be exp{dxg + icr) and exp{dyg + ic/)). Explicitly, and after 
removing some irrelevant numerical factors, the complex coordinates are 

2 -2 ■ 2 /iJ ~ J- 2ia 

z = J sm p e 

J+-J 

u;2 = i±^e2^^ (4.28) 

1 — cos P 

and the Kahler potential ends up being equal to 

/C = j_ log(j - j_) - i+ log(i+ - i) + 2j. (4.29) 

Again, a basis for the Hilbert space is given by wave functions '^\)m,n = z^vf^ . Note that, as 
in the two center case, these wave functions are, by construction, sections of a line bundle, 
C^ whose curvature is given, once more, by Q/(27r). 
To find the range ol n,m we look at the norm 

(4.30) 
This is finite if j~ < n < j^ and —n < m < n. Not surprisingly, these equations are 
identical to the original inequalities that defined the polytope, and the number of states is 
equal to the number of lattice points in the polytope; notice that this is not quite identical 
to the area of the polytope. In our case that number of points is 

M = (i+ - i_ + l)(i+ + j- + 1). (4.31) 

This connection holds quite generally for toric Kahler manifolds. As in the two center 
case fermionic contributions will correct both the state count, (4.31), and the phase space 
measure, (4.30); this will be discussed in the next section. 

The integrand in (4.30) can once more (as in the two center case) be viewed as a phase 
space density against which the supergravity solution has to be smeared in order to find 
the gravitational dual of each microstate [9] [44]. 

4.3.3 Including the Fermionic Degrees of Freedom 

From the open string point of view [14] we know that (4.31) is incorrect and that we must 
include fermionic degrees of freedom in order to account for all the BPS states (e.g. in the 
two center case). This is because, in the open string description, the centers are described 
by A/" = 4, d=l supersymmetric quiver quantum mechanics (QQM) with the position of 
each center encoded in the scalars of a vector multiplet and the latter also includes fermionic 
components which must be accounted for in any quantization procedure. 

Since we expect to see the same number of BPS states in both the open and closed 
description and since the bosonic phase spaces in both cases match exactly (and the sym- 
plectic forms agree in view the non-renormalization theorem discussed above) we may ask 
what the closed string analog of the fermions in the QQM is? 
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Consider our phase space: the coordinates, Xp, subject to the constraint (2.6), param- 
eterize the space of purely bosonic BPS solutions but, for each such solution, we may still 
be able to excite fermions if doing so is allowed by the equations of motion. If we consider 
only infinitesimal fermionic perturbations of the bosonic solutions then the former will al- 
ways appear linearly in the equations of motion, acted on by a (twisted) Dirac operator. 
Thus fermions which are zero modes of this operator may be excited without altering the 
bosonic parts of the solution (to first order). 

Determining the actual structure of these zero modes is quite non-trivial. A natural 
guess is that the bosonic coordinates of the centers must be augmented by fermionic part- 
ners (making the solution space a superspace) as is argued in [47] [48] where there is no 
potential. The fact that the bosonic coordinates are constrained by a potential complicates 
the problem in our case so we will simply posit the simplest and most natural guess and 
justify it, a posteriori, by reproducing the necessary correction to match the open string 
picture, the explicit two center and halo quantization of [14], as well as the split attractor 
conjecture [15]. 

Thus we will posit that the full solution space is actually the total space of the spin 
bundle over the Kahler phase space described in Section 4.3.2. The correct phase space 
densities are now harmonic spinors on the original phase space. This is natural from a 
mathematical point of view [49] and can be argued physically as follows. 

The space of solutions in the open string picture is spanned by letting the bosonic 
coordinates take their allowed, constant, values and setting the fermionic coordinates to 
zero (except for the center of mass degrees of freedom which we will ignore in what follows). 
We could of course try to restrict the symplectic form to this space, and then quantize, 
but this would miss the possible non-trivial topology of the fermionic vacuum. Therefore 
we will proceed in a different way as follows. 

We start with the full classical phase space of the quiver quantum mechanics including 
all the fermions. Next we are going to impose the constraints 

{Tp,H{xp)) = (4.32) 

which will restrict Xp to take values in the bosonic solution space. The constraint (4.32), 
however, is not invariant under all supersymmetries but only half of them. We can therefore 
impose (4.32) supersymmetrically as long as, at the same time, we remove half the fermions. 
The resulting system still has two supersymmetries left which one could, in principle, 
work out explicitly. If we assume that the solutions space is Kahler (which may in fact 
be a consequence of the two remaining supersymmetries), the resulting supersymmetry 
transformations will necessarily look like those of standard supersymmetric M = 2 quantum 
mechanics on a Kahler manifold. Notice that so far we have not used the symplectic 
structure for the fermions at all. 

Though it would be interesting to work this out in more detail, we finally expect 
that, after geometric quantization of the supersymmetric quantum mechanical system, the 
supersymmetric wave-functions will be £-valued spinors which, at the same time, are zero- 
modes of the corresponding Dirac operator. 
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Recall (see e.g. [50]) that on a Kiihler manifold M there is a canonical Spin'^ structure 
where the spinors take values in A'''*(A^). To define a spin structure we need to take a 
square root of the canonical bundle K = A^'^{M) and twist A°'*(A^) by that. We also 
need to remember that the coordinate part of the wave functions were sections of a line 
bundle, C. Thus altogether the spinors on the solution space are given by sections of 

£0A°'*(7W) 0^1/2^ (4.33) 

The Dirac operator is given by 

D = d + d* (4.34) 

and we have to look for zero modes of this Dirac operator. These are precisely the harmonic 
spinors on M and therefore the BPS states correspond to H^^'*>{M,C (8> K^'^). By the 
Kodaira vanishing theorem, H^^'"''{A4, C (8> K^'"^) vanishes unless n = 0. Thus, finally, the 
BPS states are given by the global holomorphic sections of £ (8> K^'"^. The only difference 
with the previous purely bosonic analysis is that the line-bundle is twisted by K^'"^. 

To find the number of BPS states we can therefore follow exactly the same analysis 
as in the bosonic case. We just have to make sure that in the inner product we use the 
norm appropriate for C K^''^. This can be accomplished by inserting an extra factor of 
{deididiK,)"^''^ in the inner product. For example, for the two-sphere, this introduces an 
extra factor of (l-|-cos 6)^^ in the integral, reducing the number of states by one compared 
to the purely bosonic analysis. This is in perfect agreement with [14]. 

For the three-center case (and more generally for toric Kahler manifolds) we find, after 
some manipulations, that 



(det9,5^/C)-V2^exp(j;|^)Jdet 



dg_ 



_cPg_ 
dxidxj 



(4.35) 



in terms of the function g given for the three-center case in (4.27). To get this relation we 
used that /C is the Legendre transform of g (see Appendix B). Evaluating this explicitly 
for the three-center case yields an extra factor 



1 /l-|-cos( 



aij) (4.36) 



j+- j y 1 - cos ( 

with 

a(j) = \/j(i+-J-)+2(i+-i)(i-j-) (4.37) 

This result indicates that, in the presence of spinors, we should take n integral and m 
half-integral with —n < m + ^ ^ n and j_ < ra < j+ — 1. Then, the total number of 
normalizable wave-functions becomes 

M={j+-j.){j++j.) (4.38) 

which does not have the unwanted shifts anymore! 
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For completeness let us provide the corrected form of the norm for a wave function, 
including fermionic corrections, 

I , |2 f -2j{3+- iy^"^ ( -2 ■ 2a(J~J- 



ij-j-V- V \3+-J 

l + cos6'\™+^ 

n «(j) J dj d cos e d(j) da. (4.39) 

1 — cos J 

Note that this is only the norm for non-scaling solutions with j_ > 0. The norm for 
wave-functions on solution spaces with a scaling point {j- = 0) is given in section 8. 

4.3.4 Comparison to the Split Attractor Flovi^ Picture 

In the previous subsections we computed the number of states corresponding to the position 
degrees of freedom of a given set of bound black hole centers. The approach we developed 
amounts essentially to calculating the appropriate symplectic volume of the solution space. 
To count the total number of BPS states of a given total charge one needs to take into 
account the fact that the different black hole centers may themselves carry internal degrees 
of freedom and that there may be many multicenter realizations of the same total charge. 
In the special case, however, when all the centers correspond to zero entropy bits without 
internal degrees of freedom the position degrees of freedom should account for all states. 
In this case it is interesting to compare the number of states obtained in our approach, 
using geometric quantization, with the number obtained by considering jumps at marginal 
stability as in [15] (see also section 2.3). 

To make this comparison we use the attractor flow conjecture which states that to 
each component of solution space there corresponds a unique attractor flow tree. Given a 
component of solution space we can calculate its symplectic volume and hence the number 
of states. Given the corresponding attractor flow tree we can calculate the degeneracy 
using the wall crossing formula of [15]. To determine which attractor tree corresponds to 
which solution space (as needed to compare the two state counts) we will have to assume 
that part of the attractor flow conjecture holds. Although this might seem to weaken 
the comparison we should point out that the attractor flow tree has no inherent meaning 
outside the context of the attractor flow conjecture thus the need to assume the latter to 
relate the former to our solutions is not surprising. Moreover, the attractor flow conjecture 
(defined in [23] [33]) is distinct from (and weaker than) the wall crossing formula (defined 
in [15]) which relies on it. 

As mentioned before (around eqn. (4.19); see also the section about the addition of 
fermions), in the two center case we get a perfect agreement between the two calculations. 
This is not so surprising because both approaches are, in fact, counting the number of 
states in an angular momentum multiplet with j = 2(^1, r2) — 5- Furthermore, there 
is no ambiguity in specifying the split attractor tree. Things become more interesting in 
the three centers case where there are now naively three attractor trees for a given set of 
centers. According to the attractor flow conjecture only one tree should correspond to any 
given solution space. It is possible to match solution spaces to attractor trees if we are 
willing to assume part of the attractor flow conjecture. 
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Let us consider the three center attractor flow tree depicted in figure 1. For the given 
charges, Fi, F2, and F3, there are, in fact, many difl'erent possible trees but, in terms of 
determining the relevant number of states, the only thing that matters is the branching 
order. In figure 1 the first branching is into charges F3 and F4 = Fi +F2 so the degeneracy 
associated with this split is |(F4,F3)| and the degeneracy of the second split is |(Fi,F2)| 
giving a total number of states 

Mree = |Fl2||(Fi3+F23)| (4.40) 

where we have adopted an abbreviated notation, F^- = (Fj,Fj) and have dropped the 
factors of r2(Fa) in (2.12) (because we are only interested in the spacetime contribution to 
the state count so we consider centers with no internal states). 

To compare this with the number of states arising from geometric quantization of the 
solution space, (4.38), we need to determine j+ and j-. As described in Appendix A, j+ 
and j_ correspond to two different collinear arrangements of the centers and, in a connected 
solution space, there can be only two such configurations. To relate this to a given attractor 
flow tree we will assume part of the attractor flow conjecture; namely, that we can tune the 
moduli to force the centers into two clusters as dictated by the tree. For the configuration 
in figure 1, for instance, this implies we can move the moduli at infinity close to the first 
wall of marginal stability (the horizon dark blue line) which will force F3 very far apart 
from Fi and F2. In this regime it is clear that the only collinear configurations are F1-F2-F3 
and F2-F1-F3; it is not possible to have F3 in between the other two charges. Since j_)- and 
j- always correspond to collinear configurations they must, up to signs, each be one of 

il = ^(ri2+Fi3+F23) (4.41) 

J2 = ^(-Fl2+Fi3 + F23) (4.42) 

j+ will correspond to the larger of ji and J2 and j_ to the smaller but, from the form of 
(4.38), we see that this will only effect M by an overall sign (the state count depends only 
on the absolute value of M) . Thus 

^f = ±(ii - J2){ji + 32) = ±Fi2 (Fi3 + F23) (4.43) 

which nicely matches (4.40). 

Of course to obtain this matching we have had to assume the attractor flow conjecture 
itself (in part) so it does not serve as an entirely independent veriflcation. Another draw- 
back of this argument is that it is not applicable to the decoupling limit described in [20] 
where the asymptotic moduli, t , are flxed to the AdS-point [20]. However, it is possible 
to circumvent this limitation by gluing an asymptotic flat region to the interior geometry. 
This can always be done by choosing the moduli in the new added region to be equal to 
the asymptotic moduli, tAdS) of the original solution. The physicality of such gluing relies 
on two important observations. The flrst one is that far away the centers behave like a 
one big black hole with charges the sum of all charges carried by the centers. The second 
important ingredient is that f AdS is equal to the attractor value of the moduli associated 
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to this big black hole. Doing so we are back to the asymptotic flat geometry where we 
have the freedom to play the asymptotic moduli once again. This is the same argument 
used in [20] [51] to generalize the existence conjecture from asymptotic flat solutions to the 
decoupled limit. 

Our result here provides another non-trivial consistency check for the conjecture that 
there is a one-to-one map between split- attractor trees and BPS solutions to N = 2 four 
dimensional supergravity [15]. By explicitly evaluating (4.31), we find it is the product 
of two contributions exactly as predicted by the split attractor flow conjecture. Using 
geometric quantization it is clear that the three-center entropy always factorizes into a 
product of splits along walls of marginal stability matching an attractor flow tree (the 
only reason we need to assume part of the attractor flow conjecture for the matching is 
to determine which particular tree). This is strong evidence in favor of the fundamental 
underpinning of the attractor flow conjecture: namely, that by tuning moduli it is always 
possible to disassemble multicenter conflgurations pairwise. 

Let us make some further remarks on the results derived here. The scaling solutions 
corresponding to A ^ have j_ = even if the centers don't align at this point. Therefore 
the connection to the wall crossing formula breaks down. The procedure of geometric 
quantization itself, however, does not seem to suffer any pathologies for these solutions. 
The curvature scales always stay small allowing us to trust the supergravity solutions. 
Thus one can see the resulting degeneracy as a good prediction. Although the symplectic 
form seems to degenerate at j = j_ this is, in fact, nothing but a coordinate artifact as 
can be seen from studying the polytope associated with scaling solutions. 

For flxed j± the Hilbert space is flnite dimensional and it is not possible localize the 
centers arbitrarily accurately. Thus the supergravity solutions can only be well approxi- 
mated in the large j± limit. In Sections 7 and 8 we will study the nature of "classical" 
states deflned in this limit. We will be interested in particular in the boundary of the 
solution spaces where classical pathologies such as inflnitely deep throats or barely bound 
centers (see [20]) moving off to inflnity may appear. We will show that quantum effects 
resolve these pathologies since there is less then one unit of phase space volume in the 
pathological regions (even for large flnite charges) so the pathologies are purely classical 
artifacts. 

5. Halo Degeneracies 

The flrst and simplest generalization beyond three centers consists of the case of so-called 
halo conflgurations [14, 15]. These are solutions consisting of a single center of charge Ti 
together with other centers that all have a charge of the form Tj = aiV^, with the aj positive 
integers^. These conflgurations are characterized by a split tree of the form {Ti, A^T^} with 
N = Ylia '^'^a, where tt-q is the number of centers of charge aF*. As all a and ria are positive 
integers every halo of total charge AF* thus corresponds to a speciflc partitioning of A^ 



*The overall sign of the a^ can of course be absorbed in F* and is thus a convention. They all have to 
have the same sign because otherwise no valid attractor flow tree exists. 
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and vice versa. It follows straightforwardly from the constraint equations (2.6) that all the 
Tj centers orbit Ti at the same distance ru = I, given by 

,_ (r.,ri) 

Note that this radius is independent of the different Oj, furthermore all the centers Fj can be 
placed arbitrarily on this sphere surrounding Ti as they don't interact among each other. 

So for a halo configuration of m orbiting centers the solution space is simply the 
product of m identical S^'s. When quantizing this system we have to take into account 
that in case some aj = aj , the corresponding centers should be treated as indistinguishable 
particles and we will have to quotient by the appropriate permutation group. As all centers 
in the Halo are independent there are two ways to proceed. Let us sketch the idea in the 
simplest case of a halo consisting of A'^ equally charged centers, i.e all a^ = 1. As we 
argued in general in the previous sections, the Hilbert space TC for a given set of charged 
centers^ is the space of sections of an appropriate line bundle on the solution space. So 
to count the degeneracies of the simple halo of A^ identical centers one could simply apply 
this recipe and construct sections of that line bundle on the solution space which, in this 
case, is {S"^)^ /S]\f. This is exactly what we will do later in this section, but actually there 
is a more standard and slightly simpler way to obtain the same result. Given that the A^ 
particles are independent and identical the Hilbert space for all the particles is nothing but 
the (anti-) symmetrized product of the one particle Hilbert space, i.e TCn = T~ii / Sm- By 
the definition of Tii this is the space of sections of Ci /Sj\f, where £i is the one-particle 
line bundle. So we can either first take the solution space corresponding to all the centers 
and construct a wave function on it, or we can take the one particle (2-center) solution 
space and construct a multiparticle wave function on it. 

In the remainder of this section we will derive formulas for the degeneracies of halos 
using both approaches and will see that the results agree. 

5.1 Degeneracy from (Anti-) Symmetrizing Wavefunctions on S' 

In this first way of counting we will consider multiparticle wavefunctions on the solution 
space of a single particle. This is very reminiscent of the approach taken in [14] . To gain 
some intuition about a Halo system consider first the simplest case of N centers, each 
carrying a charge T^,. One can think of this system as A (non-interacting) "electrons" 
surrounding a magnetic monopole of "effective charge" (ri,r^). As all the electrons are 
identical and noninteracting the wavefunction of the full system will be a symmetrized 
product of a single electron wave function. The single particle wave function is a section 
of a line-bundle C on S^, the position space of the electron. The line bundle C has Euler 
number |(ri ,r^)| (neglecting the fermionic nature of the wavefunction for a moment; see 
section 4.3.1). To generalize this to A such particles we have to tensor the bundle A^ times 
and (anti-)symmetrize due to indistinguishability. The Euler number of the resulting line 
bundle then gives the number of sections, i.e. the number of states. It is more convenient 



®Note that this Hilbert space doesn't include the degrees of freedom intrinsic to the centers, which is 
measured by their entropy. In case none of the centers has entropy then Ti. is the full Hilbert space. 
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to summarize these numbers for different A^ in terms of a generating function. Such a 
generating function for Euler numbers of the symmetric product of a hne bundle was 
given in e.g. [52]. The only difference with the discussion there is that, in our case, the 
line bundles are fermionic in nature as described in section 4.3.3. Taking this point into 
account properly, following e.g. [53], gives the following generating function 

j;d;vg^ = (l + g)Kr-r*)l (5.1) 

N 

where d^ stands for the Euler number of the A^**^ symmetric product C'^ /Sn ■ Newton's 
binomial expansion yields 

The fact that the different centers in the halo behave like fermions results in an upper 
bound of \{Ti ,r^)| for the number of such centers; this is nothing but Pauli's exclusion 
principle at work. This quantum effect is more subtle to understand from a supergravity 
perspective [14]. 

The generating function (5.3) can be generalized to include halos given by an arbitrary 
partition {ria} of N, i.e. N = X^aa?^a, where Ua is the number of centers carrying the 
same charge Ta = aT^,. It is not hard to see that the generating function including such 
arbitrary partitions is given by 

E^^'?^=n(i+^')"^'"'*^ (5-3) 

N fc>0 

where D^^ is the degeneracy of all halos with total charge NT^. The degeneracies can by 
found by expanding the product: 

where the sum is over all possible partitions P = {ua} of A*". This agrees with the fact that 
the total degeneracy of a given partition is the product of the degeneracies of each group of 
identical terms in the partition. The degeneracies of such a group of identical halo charges 
was exactly what we calculated in (5.1) to be a binomial coefficient. 

The formula (5.3) is the similar to (5.6) in [15] but there are also some obvious dif- 
ferences. The first one comes from the fact that [15] is calculating an index while we are 
counting states without relative signs. The second difference is that we are neglecting the 
degeneracies associated to internal degrees of freedom of each individual center. 

5.2 Degeneracy from Wavefunctions on (S^)"/5n 

We can now reproduce the above results in a more straightforward way using techniques we 
developed in earlier sections. Here we will directly quantize the solution space correspond- 
ing to a given halo. As discussed before, the solution space for halos will be the product 
of various two center spaces, i.e. spheres, quotiented out by the appropriate permutation 
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group, due to the identification of identical centers. An elegant way to describe this space, 
that can also be generalized to more involved configurations of the next section, is by using 
the technology of toric orbifolds that was developed in [54, 55] (see Appendix B). 

Let us start again with the simplest halo configuration, that of N centers of charge T^ 
orbiting a charge Ti. As discussed in the beginning of this section the classical solution 
space for this configuration is the product of N S^'s. The polytope describing each of those 
is simply an interval, see e.g. [46], 

-j < Xi< j 

where j = ||(ri,r^)| is the size of the two-center angular momentum and Xi = j cos 9i. 
The polytope associated to the product (S^)^ is then an A^-dimensional hypercube. 

Because we are dealing with N identical centers we need to quotient the total space by 
the permutation group S^. The resulting space is a toric orbifold [54]. Such toric orbifolds 
are uniquely described by labelled polytopes. These are polytopes where each facet, i.e. 
(d — l)-dimensional face, comes with a positive integer m, where the points on the facet 
are Z^ singularities. In our case of S^ /S^ the polytope is given by the convex hull of the 
following facets with their corresponding labels: 

• two facets with label 1 given by 

i) the facet —j = xi < X2 < ••.• < xj\[ < j. 
ii) the facet — j < xi < X2 < ■■■■ < xn = j- 

• {N — 1) facets with label 2 given by 

iii) the facets —j < xi < .... < xi = Xj+i < .... < xn < 3 (one for each i = 
l,...,Ar-l) 

This polytope is the quotient of the original hypercube by the orbifold group Sn and the 
labels can be easily understood by considering the fixed points of the group action. 

Given a labeled polytope there are canonical complex coordinates on the corresponding 
toric orbifold [55]. As in the non-orbifold case they are constructed from a single function 
g which, in our specific case, is given by 

1 1 ^"^ 

g{x) = -{xi +j)ln(xi + j) + -(j -a;n)ln(j - Xn) + ^ (xj+i -Xi)ln2(xi+i - Xi) . (5.5) 

j=i 

The complex coordinates are then given in terms of this function as exp(9j(7(x) + i^j) 
leading, in the above case, to 

4 ~ - , 7"'"' . <?'*' (5.6) 

J (cos6'2 — coso'i)"^ 

^, ^ ^cos^^-cos^y ^2.,. ^ ^ = 2,...,iV-l (5.7) 



cos Wj+1 — COS t/j 
(cos On — COS 6 IS 

j 1 — COS On 



Zl ~ i ^^^-Qn-^O^Q^N^^f ^2.,. (58) 
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These are the complex coordinates for which the symplectic form (4.10) is a Kahler form. 
A basis for the Hilbert space of holomorphic sections of the hne bundle with ci(£) = 
rj/27r, is then formed by the wave functions "^m = Y\i=izT'- '^^^ range of m = (mj) 
is constrained by requiring that the tpm are normalizable. Due to the fermionic na- 
ture of the wavefunctions the appropriate integration measure to calculate the norm is 
e~^(det 9icfc/C)~^/^ (as in section 4.3.3). The Kahler potential can be calculated from g 
(see [55]) and is, in our case, 



^ ~ — j ln[(l + cos^i)(l — cos^at)] + j{cos6n 
With some manipulation we determine that 

dg 



cos^i) 



(5.9) 



(det5.5^/C)-V2~exp Q:^ detV^ 



d'g 



dxidxj 



and furthermore use that 

N 

exp 



E 



dg 
dxi 



1 + cos 6i 
1 — cos On 



(5.10) 



det 



d'g 
dxi dxj 



1 



n-l 



n 



1 



[1 H (cos ^7V — cos 6*1)] — /111 

2 1 — costal 1 — COSt'Ar -'■-'■ COS yj+l — COS ( 



Putting everything together we find the following norm for the wavefunctions 



\A, 



f e-^(yN-yi)[l + l(cos0;v - COS0i)]^/2(i ^ ^^y+mi^^ _ ^^y-m^-i 



(y2-yi)'('"^-'"^^-'/' {yn-yN-i, 

Normalizability then implies the following domain for m: 



2(m„-mjv_i)-l/2 



(5.11) 



-J < 'iT-i < m2 < .... < niN < j 



(5.12) 



This can be interpreted as picking a representative of a totally antisymmetric A^-tensor, 
which is the expected behavior of N identical fermions. The degeneracy is thus nothing 
other than the binomial coefficient, which means that the generating function is given by 



^dNQ 



N 



:i + 



il(ri,r.>| 



(5.13) 



N 



where we have used the fact that the angular momentum j is proportional to the intersection 
product of the charges. So we find the same result as in the previous approach (see (5.1)). 
We can now also extend this construction to more general halos. One can reduce 
the more general case given by an arbitrary partition of N, to products of our previous 
simpler example by grouping all equal terms in the partition. Writing partitions of N as 
N = ^^ aua, the solution space is then of the form Yl^i^^)"'"' /^^a- The derivation of the 
complex coordinates and the normalizability constraint then proceeds analogously to the 
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previous example. Using the indices m" = (m^), i = 1 . . .Ua to label the wavefunctions 
V^m" = riiC'^j)'"' °^^ finds that they are normalizable iff 

-aj <mi <m2< .... < m^ < aj (5.14) 

The degeneracies of all the halos with halo charge NT^ can then be captured using the 
following generating function: 

E^^'^''=n(i+^')"^'"'*^' (5-15) 

AT fc>0 

which coincides with eqn. (5.3) again showing the equivalence of the two derivations. 

Using this generating function we can derive the large N behaviour of the degeneracies 
Dj\f. Let us define 



I -lyt—i p— /J" 



N n=l 



where in the last equality we used expression (5.15) and some rather straightforward ma- 
nipulations of the resulting series. The physical interpretation of F and /? is that of the free 
energy and inverse temperature. As is familiar from thermodynamic systems we can find 
the large energy, or in this case large N, degeneracies by approximating the free energy at 
large temperatures and then Legendre transform to find the corresponding entropy. In the 
large temperature limit /3 <C 1 the free energy is approximately 

F(/3)«^Kri,r.)|c(3)^, (5.17) 

with ^(3) ~ 1.202 the Riemann zeta function evaluated in 3. From the free energy we can 
calculate the average "energy" N in terms of /? by the standard formula 



n 77 O 1 

5^ = 2Kr..r.)IC(3)^ 



« = -»:j = :;Kri,r.)IC{3)^^ (5A8) 



Note that indeed the large temperature regime corresponds to that of large N. Using the 
relation between free energy and entropy we can then find an expression for the large A^ 
degeneracies: 

3 ,,^ ^ .,._,.2^^ 



SiN) = log DN = m-F=i-\{TuT^)\ as) N'\ , (5.19) 

where we used the relation (5.18) between N and /9. The conclusion is thus that for a 
halo of total charge NT^, there are of the order A^" states when A ^ 1. Note that this 
is slightly larger than the growth of A"*^" that comes from the partitioning of A, due to 
the presence of additional angular momentum degrees of freedom that also grow with A. 
As the angular momentum is due to the electromagnetic coupling of the halo particles this 
"entropy enhancement" is in some ways similar to that observed in e.g. [56] and [57]. 
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6. Dipole Halos 

So far we have restricted our analysis to generic solutions with two or three centers and a 
first generalization to an arbitrary number of centers in the specific case of a single charge 
surrounded by a halo of identical centers. The next further natural generalization is to 
bind the halo to an additional center. More precisely we can embed the halo split tree 
{Fi, A^r^} in a larger tree of the form^'' {r2, {Ti, A^F^}}. If the moduli are very near the 
wall of marginal stability the additional center T2 will be very far apart from the halo and 
the degeneracies will simply factorize [15], i.e. Dj\f = ((r2,ri) + N{r2,Ti,))D^^°. Note 
that to derive this degeneracy we use the attractor flow conjecture which only works for 
non-scaling solutions so for scaling solutions we will have to resort to other methods. For 
more arbitrary values of the moduli the center F2 will be closer to the halo and deform 
it. In certain cases it can even deform so much that the topology of the split changes to 
{Fi, {F2, iVF^}}. Such a change can happen when crossing a wall of threshold stability 
[20] [15] and in that case the number of states does not change (even though the topology 
of the tree changes). 

In this section we will reproduce the degeneracy as inferred from the simple flow tree 
arguments above by our quantization methods. We will do this in the specific example of a 
D6-D6-D0 system, which seems closely related to the microstates of the D4-D0 black hole 
[58]. There exist such configurations with a purely fluxed D6-D6 pair and an arbitrary 
number of anti^^-DO's. Depending on the sign of the B-field these DO's bind to the D6 
or anti-D6 respectively. When we take the B-field to be zero at infinity the system is at 
threshold and the DO's are free to move in the equidistant plane between the D6 and anti- 
D6. This system and its behaviour under variations of the asymptotic moduli was studied 
in detail in appendix B of [20] . 

As the number of states will be independent of the asymptotic moduli (as long as 
we don't cross a wall of marginal stability) we are free to choose them such that the 
solution space has its simplest form. The symplectic form on solution space is most easily 
calculated at threshold which, in our example, corresponds to i?|oo = [20]. For a set of 
DO's with charges Fa = {0, 0, 0, —Qa} with all the qa positive and J2a ?« ~ ^' bound to 
a D6, Fg = (1, 1) ^) fg)) and D6, Fg = (—1, |, — ^, fg), the constraint equations take the 
following form at threshold: 

_^ + ^ = (6.1) 

-— + E— = -/5 (6.2) 

3 

Here / = — (F6,Fg) = ^ is given in terms of the total D4-charge p of the system and 
/? = {TQ,h) with 1,(3 > 0. From the first line we indeed see the DO's lie in the plane 
equidistant from the D6 and D6, as we are at threshold, and so we can simply write 

^"in this notation figure 1 would be {T3, {Fi, r2}}. 

^^In our conventions it is anti-DO's that bind to D4 branes. We will however often just refer to them as 
DO's. 
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An explicit expression for the symplectic form (4.10) can be obtained using the fol- 
lowing coordinate system. We define an orthonormal frame {ii, v, w) fixed to the D6-D6 
pair, such that the D6-D6 lie along the w axis and with the DO's lying in the u-v plane. 
Rotations of the system can then be interpreted as rotations of the {u,v,w) frame with 
respect to a fixed {x,y,z) frame. We will parameterize these overall rotations in the stan- 
dard fashion by two angles, {9, (p). We can furthermore specify the location of the a'th DO 
with respect to D6-D6 pair by two additional angles, {Oa,(j)a)- The first angle, 9a, is the 
one between Xgg and X6a, while (j)a is a polar angle in the u-v plane (see figure 3). Our 
2N + 2 independent coordinates on solution space are thus {9, 4>,9i,(j)i, . . . , 9^-, (J^n}- 




Figure 3: The coordinate system used to derive the D6-D6-N DO symplectic form. The coordinates 
{9, (j)) define the orientation of the {u,v,w) axis with respect to the fixed, reference, {x,y,z) axis. 
The D6-D6 lie along the ii axis (with the origin between them) and the DO's lie on the ii-v plane 
at an angle (j)a from the il-axis. The radial position of each DO in the ii-v plane is encoded in the 
angle 9a (between Xgg and a?6a)- 

The standard Euclidean coordinates of the centers are then given in terms of the 
angular coordinates by 



xe 



J_ 



w 



j tan^a 
Xa = "5 (cos aau + sm aav) 



u = COS (p X — sm (j) y (6.3) 

V = sin (/) cos 9 X + cos (j) cos 9y — sin 9 z (6.4) 
w = sin (j) sin 9 x + cos (j) sin 9y + cos 9 z (6.5) 



The angular momentum, j{9a,4'a), is a function of the other coordinates rather than an 
independent coordinate (when A^ = 1 it can be traded for 9i as demonstrated in the general 
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three center case), and is given by 

J = ^-^Qacosea. (6.6) 

a 

Using this exphcit coordinatization it is straight-forward, though tedious, to evaluate the 
symplectic form (4.10). The result is relatively simple: 



4 



2j cos 9 d(j) + 2 2^ Qa cos 9a da a 



(6.7) 



with d denoting the exterior derivative. Note that, as is manifest from our angular co- 
ordinatization, we are still in the toric setting with each additional center introducing an 
additional U{\) coordinate. Note also that for a single DO this reduces to (4.25) when 9i 
is traded for j. 

In case N > 1/2 it is possible to combine a sufficient number of centers and form a 
scaling throat. When there is more than one DO center this spoils the simplicity of the 
polytope and the toric manifold becomes singular. This is one of the reasons for which we 
will restrict our analysis in this paper to the non-scaling case 1/2 > N , or, as in section 
8, restrict to three center scaling solutions. We hope, however, to return to more general 
scaling solutions in the near future. 

6.1 Degeneracy using Toric Techniques 

In the following we are going to use the same toric techniques as in the previous sections 
to calculate the degeneracy associated to the non-scaling D6-D6-D0 states, i.e. those for 
which A^ < 1/2. From the associated symplectic form, eqn. (6.7), one easily reads the toric 
coordinates to be 

y = j COSO , Xa = Qa COS 9a > (6.8) 

Because some of the centers can be identical we need to orbifold our polytope by the appro- 
priate symmetric group. Before treating the configuration given by a generic partitioning 
{qa} of N i.e. N = ^^ naQa, let us focus on the simple case of n centers carrying the same 
DO charge, —q, so that N = qn. In this case the polytope without identifications is given 

by 




(6.9) 

where we have used the explicit form of the angular momentum j from eqn. (6.6) in the 
first inequality. We can then orbifold this polytope and introduce labels to the facets of 
the quotiented polytope. In the current example this amounts to 

• four facets with label 1 given by 

i) the facet xi = 0, 
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ii) the facet Xa = q, 
iii) the facet y = {-Y^a^a, 
iv) the facet y = -5 + Xla ^a, 

• (n — 1) facets with label 2 given by 

v) the facets Xa+i — Xa = 0, for a = 1, ..., n — 1. 

Given this labelled polytope we can then again construct the corresponding complex 
coordinates through the function 



^ a^^ a^^ a^^ a 

1 1 """^ 

+ -xilnxi + -{q - Xn)\'a.{q - x„) + ^(x^+i -Xa)ln2(xa+i - Xa) 

a=l 

Explicitly they can be calculated by Zi = exp (dig + i(j)i), leading to 

° //2-y-Ea^a' 

z? ^ I ^ e^'"' 

4 ~ 7777^ F^ ^Wtt^ F^ r ^'':~''"'^^' e^-" (6.10) 



(//2 + y - Ea ^a) (//2 -y-Ea ^a) (q " X„) 



.2 1 / ^ ^llzi '\ ' ^2ia. 



2 



<~777;7T Y^ ^777o v^ V ^ ^ ^ ^ ^ = 2,...,n-l 

[I/2 + y-X,^Xa){I/2-y-l^^Xa) \Xi+i -XiJ 

The next step is to construct the Kahler potential, which turns out to be equal to 

1^= -^'^^U+y-T.a^ajU-y-T.a^a] 

-2J2a^a + Xn- xi- qln{q- Xn) (6.11) 

A basis for the Hilbert space is then given by normalizable functions ipm = Y\7=o -^I"' where 
the norm is given by 

As in our previous examples the determinant can be computed by the identity 



'detdidjlC ~ exp[9jg'] ^ydet^i^jg . 
The first term is simply Y\- \zi\. The second term is more involved but can be shown to be 



1 / T \ -1 "-1 



detdidjQ r^ i- + y-'^Xaj [2~y ~'^^'') ^^^i) ^(Q-Xn) ^ f| (^a+i - x^) ^ 

^ a ^ ^ a ^ a=\ 
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Using these expressions to analyze the normahzabihty constraint (6.1) then leads to the 
following constraint on the possible exponents m = {rrii): 

< mi < 1712 < .... < rUn < q , 

-i^-Za {ma + l))<mo + l<{^-Za K + D) (6-12) 

where n is the number of DO-centers carrying charge q. The total number of normalizable 
states is thus 

(/-n-l)/2 

dn,,= Y. b"miQW-n)-2m] (6.13) 

m=n(n— 1)/2 

where the coefficient b'^{q) is the number of ways to write m as a sum of n strictly ordered 
positive integers all smaller than q. 

Let us now generalize the simple example of n equally charged DO-centers to an arbi- 
trary partition of N. Label the different groups of equally charged centers by an index a 
and the charge of individual centers in this group by qa (i.e. (7^ 7^ % iff a / 6). With Ua we 
then denote the number of centers with charge qa so that the total DO-charge N carried 
by the DO centers is given by 

A^ = ^ naqa 

a 

Labeling the centers in a given group a by an additional index i = 1 , • • • , ?!„ we can simply 
generalize the conditions (6.12) by applying them to each group of equally charged centers 
separately. The conditions on the powers m"' = (?n-f ) then become 



< mf < m| < .... < m^^ < qa , 



l-i 

2 



Ea,^« + 5) <mo + l< 



7-1 



Ea,i « + 2) 



(6.14) 



A first step towards counting all possible states with total charge A^ in DO-centers is to count 
the degeneracy for a fixed partitioning of N. The number of solutions to the constraints 
(6.14) is given by 



d 



{na,qa} 



E 



all allowed rra? 



/-2E(<+^) 



a \ "■ / all allowed mf i,a ^ ^ 



(6.15) 



We can calculate the sum of the last terms by introducing the quantities 



l'} = qa-l- m" 



and noting that then 



< /^ < /^ < 



< ^n„ < la ; 



(6.16) 



Z^all allowed mf Z^a,i ("^i + 2) ~ ^ Ha 



ria 



X/all allowed If z2i,a {^i + 2) V^-^''') 
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where we used that Ylia'^^a^a = N. As /" and m" satisfy the same conditions, eqn. (6.17) 
simply implies that 



all possible m" i all possible '" * 



Using this last equality we see that degeneracy (6.15) is given by 




d 



{na,qa} 



,.-.)n(r>u-E»«*)n(: 



(6.18) 



So we indeed find back the result we derived using attractor tree arguments: the degeneracy 
is that of the corresponding halo, multiplied by (ri,r2) + A^(ri,r^) = I — N. 

Counting all the different degeneracies for all possible partitions of a given total halo 
charge NT-^, gives rise to the following generating function: 



Z{q)=Y,DNq'' 



N 



{na}{qa} « a 




ria Qa 



iI-Qd,)ll 



Qa 



E 



^^ I „naqa 



n,, 



[I - qd,)\{{l + q^)K 



(6.19) 



Using this generating function we can again estimate the large A^ growth of Djv, but, 
as the only difference from the (non-dipole) halo case is a pre- factor of order {I — N), we see 
that the growth of log -D at ~ N'^'^ also holds in this case (modulo logarithmic corrections). 
Note that, although we find an exponential number of states, the N'^'^ scaling of the entropy 
is far less than the expected N (at A ~ /) for a black hole with charges //4 < N < 1/2 (for 
N < 1/4 the state is polar and no single center black hole exists with this total charge [15]). 
It would, however, be extremely interesting to do a similar counting for scaling solutions 
{N > 1/2) which do admit single center black hole realizations and compare this to the 
black hole entropy. We hope to come back to this in the near future. 

7. Quantum Structure of Solutions I: Centers Near Infinity 

Having determined the structure of the quantum states associated with the two and three 
center solutions we can now investigate some potential problems with the classical solutions 
that we expect to be resolved by quantum effects. For instance, although we have already 
argued that the apparent non-compactness of the three center phase space for some choices 
of charges and moduli is not a problem since the symplectic volume of the phase space is 
finite we would still like to see that centers cannot move off to infinity once quantum effects 
are taken into account. Likewise, as has already been observed in [35] and [15], scaling 
solutions can develop an infinitely deep throat classically but we expect the quantization 
of phase space to cap this throat off at some finite value. 
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To investigate these issues we would like to consider the expectation value of the 
harmonic functions (2.3) defining the solutions. Of course this will depend on the particular 
state we are considering and there are, in general, many possible states one can construct 
so making any general statement is quite difficult. We do not, however, need detailed 
properties of {H{r)), only its behaviour in various asymptotic limits. In this section we 
are going to discuss the non-compact case leaving the scaling solutions to the next section. 
Our main concern here will be the r dependence at large r, which we will be able to extract 
for a general pure state. 

Let us consider the two independent constraint equations (3.2) for the three center 
case once more 



a 
u 


b 

V 


ri2 

ri2 




= {h,T,) 


b 

V 


c 

w 


_r3i 

^-31 


r23 


= -{h,T3 



--: a 



--■■P 



(7.1) 

(7.2) 



where we have (re)introduced a hopefully obvious short-hand notation. We would like to 
see when one center can move off to infinity. We will assume that a, b and c are non-zero. 
In order to satisfy the triangle inequalities while taking at least one center to infinity we 
must have either two or three of the distances u, v, and w become infinite. Let us, for 
definiteness, try to set u and v to infinity which corresponds to centers 2 and 3 staying 
a finite distance apart while center 1 moves off to infinity. From the constraint equation 
its clear that this can only be done if a = 0. We can also consider the case when all 
three centers move infinitely far apart, which can only occur if /3 = as well, but this is 
a somewhat trivial case as the angular momentum will then, by (2.8), vanish, as will the 
symplectic form, so the space cannot be quantized without adding additional degrees of 
freedom (momenta). 





Figure 4: (Left) Possible locations for a?i are given above. The centers X3 and .T2 sit at (0, 0) and 
(0, —2) in the (x, z)-plane respectively and the different orbits represent different values of a/b. hi 
the figure above a/b < 1 and as a/b — > 1 the orbit increases until it becomes the line z = — 1 when 
a = b (not depicted). (Right) Normalized probability densities \ipm.n{j)\ plotted as a function of 
j alone (neglecting angular dependence) for some values of n with J^ < n < J^. In this example 
J+ = 18 and J_ = 6. 

Thus we restrict to the case a = 0, /? 7^ which gives 
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- = - = ? (7.3) 

rsi V b 

When translated into coordinates this equation defines an ehipse for the possible locations 
of center 1 but it is easy to see from this that u and v must be bounded unless a = 6, in 
which case the orbit is not an ellipse but rather center 1 must lie on the plane between 
centers 2 and 3. Hence all unbound three center solutions will be of this form (where center 
1 is fixed on the plane normal to the axis defined by the other two centers). Orbits for 
various values of a/b are depicted in figure 4. Note that this proves the claim of [20] that 
for three centers 'non-compactness' can only appear when the asymptotic moduli are at 
threshold stability, as a = 0, a = 6 is exactly the definition of threshold stability [20] in 
different notation^^. 

Recall that J^pi^^^p) = so if we define 7 = {h,r2) then a = implies /? = 7. Using 
(2.8) we find that 

J = -r23Z = —z (7.4) 

with the positive z axis defined by X3 — X2- 
Let us consider the other equation 

b c 

---=/3 = l (7.5) 

V w 

where we rescale the coords so /3 = l}^ Note that this forces c < if we want to allow 
u — > 00. 6 can be positive or negative but we will assume 6 < for concreteness (this will 
not alter the analysis). 

Using the fact that r23/2 < ri3 < cxd we find 

J+ = y i-=Min(^-|6|,o) (7.6) 

with j+ reached at ri3 — > cxo and j- at ri3 = |c|/2 — |6| (unless this is less than zero^^). 
We can now consider the expectation value of 

H{r) = -K- (7.7) 

\r — ri\ 

The asymptotic behaviour of this function is particularly important in the decoupling limit 
considered in [20] since if ri ^ co the classical solutions will no longer be asymptotically 
AdSs. Thus we would like to check if there are wave- functions for which {H{r)) does 
not decay as r~^. If such states exist they would spoil they asymptotics of our solutions, 



^^More precisely: a = 6 <;=> (Fi, Ta + Ta) = and a = <^4> lm{ZiZ2+3) = 0. 

^^We can permute the centers to force a positive sign since (3 appears in one of the three contraint 
equations and —f3 on the other and we are free to use any two of the three. 

^^Since a = b, b — c/2 corresponds exactly to the beginning of the scaling regime and is, in fact, nothing 
more than A'' = J/2 in the specific example of D6D6D0. Although we work here with j_ > it should not 
matter much since that large n states we consider in this section have little support in the small j regime. 
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particularly in a decoupled AdSs limit as in [20], and it would be hard to interpret them 
physically. 

We are interested in studying wavefunctions localized near infinity so we could restrict 
our attention to the states with the largest angular momentum, n = j^ — 1, but it turns 
out to be tractable to study more general pure states cj) = J2n m '^n,mipn,m (though we will 
find contributions from n < j^ — 1 are more strongly suppressed near infinity as suggested 
by figure 4, right). For this state the expectation value is given by 



POO 

n'm\n,m i 



n' ,m' ,n,m 



{ri) dr I d COS 9 d(j) da (7.^ 



h{r,ri) = — (7.9) 

(r^ + ri — Irri cos a(y, 0, o"j j ' 



2N+1 



(7.10) 



V(Aj+)(j+ - j„) + 2j+(l - A)(Ai+ - j„) ^^^'^l ^^^ (7.11) 

A(.,) = ^ (7.12) 

Here N = n + n' so j_ < N /2 < j+ ~ 1 ^^d n_ = n_3 because x^, is at the origin. The 
function fn',m',n,m = Ki' ,m' '^ri,m i'ni ^rn' '^n,m dj / dr which we re-write in terms of ri using 
j = A(ri) J+.^^ Since we are only interested in the large r behaviour of this function we 
can integrate out the angular dependence and also neglect constant factors (both of which 
have been absorbed into the function g). 

The integral splits into two parts given by r > ri and r < ri. Since A(ri) r^ l+b/ri + . . . 
for large ri, in the second region we see that /(ri) ~ r-^ so, after performing the 

angular integrals, we are left with 

f: C,,,/ r r-^-'-'-^'--"^/'^ (7.13) 

l,k=0 •^'" 

The expansion in k comes from expanding A(ri) in powers of r^ . Clearly the ri > r region 
only contributes negative powers of r to {H{r)). 

In the region ri < r we cannot expand A(ri) since ri may not be larger than \b\ but 
we can split this integral once more into two regions: r_ < ri < f and f < ri < r for some 
f S> |&|. In the first region the integration domain is r-in dependent so the r dependence 
is simply r^'~^. In the second region we can repeat the analysis for the ri > r integral, 
expanding A in rj~ , and we find a similar greater than rj~ fall-off. 

As was mentioned above, such configurations (with a = and a/b = 1) lie on walls of 
Threshold Stability as discussed in [20, Appendix B]. In fact, our computation nicely agrees 



^Note that we have absorbed a factor of (ri — 6) ^ from the jacobian dj/dri into our definition of f{ri). 



40 



with the fact that, as follows from the wall crossing formula [15], no states actually decay 
when crossing such a wall of threshold stability. More loosely speaking, our calculation 
roughly excludes the possibility of "states running off to infinity". This is important e.g. 
in the consistency of the decoupling limit discussed in [20] where the limit itself forces the 
moduli to a wall of threshold stability for many charge configurations and we do not want 
this to spoil the asymptotics of the solution. It is also more important in a more general 
context as unbounded centers do not admit an easy physical interpretation. 



8. Quantum Structure of Solutions II: Scaling Solutions 

We have seen in the previous section that the naive run away behaviour of some solutions 
is, in fact, a classical artifact which vanishes once quantum effects are correctly accounted 
for. The other potentially paradoxical region in solution space is near a scaling point, where 
an infinitely deep throat develops in spacetime [35, 15, 17] (see also section 3.2.1). This 
throat is particularly enigmatic in the context of AdS/CFT. An infinite throat suggests 
that the bulk excitations localized deep in the throat will give rise to a continuum of states 
in the dual CFT. 

To make this connection with the dual N = (0,4) CFT we uplift the 4-dimensional 
solutions to 5-dimensions and take the decoupling limit described in [20]. The quantization 
procedure described so far will carry on mutatis mutandis to the uplifted solutions because 
they have exactly the same solution space. We will start with a general discussion but then 
specialize to a working example given by a three center D6D6D0 scaling solution where 
the charge of the DO center, N^ satisfies A^ > (r6,rg)/2. We begin with some details on 
the structure of scaling three-center solution spaces since, as alluded to in section 4.3.2, 
there are some subtle differences in the geometry of scaling and non-scaling solution spaces. 
After constructing the appropriate wave functions from the Kahler geometry we will use 
them to estimate the depth of the throat. We will argue that these throats get capped at 
some scale e and that this also sets the mass gap for the CFT. We will perform an estimate 
indicating that e ~ A^/(r6,rg) and argue that the corresponding mass gap in the CFT has 
the familiar 1/c behaviour when e is of order 1. 

8.1 Quantizing the Three Center Scaling Solutions 

As was done in the non-scaling case we must first construct the appropriate polytope for 
these solutions (see figure 2). The only property that differentiates these solution spaces is 
that j- = (this is the scaling point). As a result, the associated polytope differs slightly 
from the non-scaling one; for instance, the first inequality in (4.26) is redundant. This 
may seem to be a small modification but it actually changes the topology of the solution 
space, taking the limit from non-scaling to scaling corresponds to a blow down and a non 
contractible 5"^ vanishes. Furthermore as we will discuss later, the probability densities at 
the boundary of solution space, j = j_, will have a very different behavior in the j_ = 
case (figure 5) than in the j_ / case (see figure 4). 
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Using the coordinates x = j and y = jcosO as in section 4.3.2, the scahng solution's 
polytope is defined by 

j+-x>0, x + y>0, x-y>0 (8.1) 

The construction of the complex coordinates is achieved through the function g (see ap- 
pendix B). We will only need an expression for their norm squared, which is given by 

, ,9 7^sin^6' , ,9 l+cos^ , , 

^if = - ^, 22' = . ^. 8.2 

J+ - J 1 - cos fc* 

To get the wave functions one requires that they have a finite norm using the measure e~^, 
modified by fermionic corrections as discussed in section 4.3.3. /C, as usual, stands for the 
Kahler potential. One can calculate that 

K=i-i+log(j+-i), (8.3) 



de.(8.8;r) = w'l±^^^5^3^ (8.4) 

■' ' V 1 — COS & J+ — J 

Putting everything together, we end up with the following norm 

(1 + COS ^)"+('^+l/2) (1 _ cos e)n-[^+l/2) ^- ^^^^ q (35) 

Requiring that the norm is finite imposes the following restrictions 

< n < j+ - 1 , -n < m + 1/2 < n (8.6) 

So the number of states is given by 

Unfortunately, we cannot compare this prediction to wall-crossing because it is not clear 
how to treat scaling solutions within the framework of the attractor flow conjecture [15]. 
On the other hand this proves the usefulness of the tools we developed in this paper as they 
provide the only known way to compute the number of BPS states for scaling solutions^^. 
Another important property that is worth mentioning is that the probability density, 
given by the integrand of (8.5), vanishes at j = 0. This suggests that, although classically 
the coordinate locations of the centers can be arbitrarily close together, quantum mechani- 
cally this is not true any more. The probability that centers sit on top of each other is zero 



^^While the degeneracies computed for scaling solutions cannot be checked against the attractor flow 
conjecture there is, nonetheless, a tantalizing connection with a particular open string computation. If we 
consider a configuration with a pure fluxed D6 and D6 center with charges {±l,p/2, ±p^/8,p'^/48} and a 
single DO-center as in [58] (see also section 6 where the non-scaling version of these configurations were 
considered) then the number of states, A/", is proportional to p®. Additional DO centers will increase the 
state count by a factor of p^ each. This is very reminiscent of [56, 58] where orbiting DO's puff up into D2's 
which can each occupy one of ~ p'^ Landua-levels in the Calabi-Yau. We hope to make this connection 
more precise in the near future and see if it can lead to a derivation of black hole entropy in terms of zero 
entropy supergravity solutions. 
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which imphes that there is a minimum non-vanishing expected inter-center distance. Since 
the depth of the throat is related to the coordinate distance between the centers it fohows 
that the throat will be capped off once quantum effects are taken into account. In the 
following section we will study this phenomena quantitatively and make some predictions 
for the depth of the throat and the corresponding mass gap in the CFT. 

8.2 Macroscopic Quantum Effects 

Before we analyse a specific type of scaling solution in some more quantitative detail let 
us point out some general features of any three center solution space with a scaling point. 
As long as we are interested in spherically symmetric quantities, e.g. sizes and distances, 
we can neglect the angular part of (8.5) as it will drop out after normalization. For such 
questions we can effectively use wave functions on j-space, which depend only on the 
quantum number n and the value of j+ . The corresponding probability densities are 



{n,j+;j\n,j+;j 
These are plotted for n = 0, 5, 10, 15, 20 and j 



Jo 



e-if^+\j 



^-jj2n+l^j^ 



21 in figure 5 (left). 



(8.7) 





Figure 5: (Left) Plot of the probability densities (8.7) for j+ = 21. The blue, orange, red, purple 
and green curves correspond to respectively n = 0, 5, 10, 15, 20. (Right) Plot of the probability 
density corresponding to the lowest state n = for the values j+ = 5, 50 (blue) and j_|- = oo (green) 
where the curve at j+ = 50 is already barely distinguishable from the limiting curve j+ = cx3 (see 
eqn. (8.9)). Note that the probability distribution vanishes on j = 0. 

The state that interests us the most is the n = state as this has the greatest support 
near the classical scaling point, j = 0. Note that when there is no scaling point, i.e. j- ^ 0, 
the lowest state, with n = j-, peaks on j_ (see figure 4). In the scaling case the lowest 
state actually has zero support on j_ = (see figure 5, right). This seems to indicate that 
the scaling solutions, supergravity solutions where all centers coincide in coordinate space 
and j = 0, are not well defined classical solutions as they correspond to a point in the 
phase space where no wave function has finite support. To get an idea of how well this 
classical point can be approximated by a quantum expectation value we calculate {j) in 
the lowest state, |0, j+;j)- 

For a scaling solution space characterized by j+ this expectation value is given by 



ij) 



3+ 



■J+ 



{0,j+;j\j\0,j+;j)dj 
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In general this expression is not analytically tractable. We are, however, particularly in- 
terested in the supergravity regime which coincides with j_|_ — > oo. In this limit the expres- 
sion can be simplified considerably by using the well know expression for the exponential, 

limj_^^oo(l - ^y+ = e~^+, giving 

.lim (0,i+;i|0,i+;j)=4e-2^j, (8.9) 

this is the green curve plotted in figure 5 (right). Using this limiting behaviour it is 
straightforward to calculate that 

(j)oo = 1 . (8.10) 

In other words, even in the lowest state the expected value of j is one quantum, i.e h. 
Moreover, because the depth of the throat grows very rapidly in the region j ~ many 
macroscopically different configurations sit within the range < j < 1 so, even though 
J ~ 1 is only one plank unit away from the scaling point the corresponding geometry 
(expectation value of the metric) is very different. 

As further evidence for the formation of a cap we will also compute the behaviour of the 
harmonic functions appearing in the metric in this same state. To make this computation 
simpler we will work with the system introduced in section 7 (a = and a = b but now we 
take |5| > |c|/2) and we will also use the notation used in there. We will first be interested 
in determining the r-dependence of the (expectation values of the) functions 

hi{r) = ^^ h2{r) = — ^ (8.11) 

|r — ri| \r — r2\ 

where ri and r2 have the same meaning as in section 7 (rs = — r2 so we need only 
work with one of them). Note that the functions hi{r) appear directly in the harmonics 
{H^{r), H (r), Ha{t)., HQ{r)} but which harmonics they appear in depends on the specific 
form of Ti (which we do not fix at this point). 

To compute the r-dependence of hi{r) we proceed very much as in section 7, and indeed 
the computation is mostly analogous, 

h.(rr)-i S£oqo(^''^''^)4 "^''' (8 12) 

i 

r^O,j+;i|/i,(r)|0,j+;j)dj = i V[r-'-i f C^i)r\ f{n) dn+ 
JO ^iiJo 

r' f C[i^ r7^-' f{n) dn + r' H C|,) r"'"! f{n) dr} (8.13) 



f{n) = j{n) e-2^('^»)V2J+-j-(n)^ (8.14) 

i(n) ~ -^\j+ (8.15) 

ri — 

j{r2) ~ r2 (8.16) 
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Here we are only interested in the regime r <C |6| and we define f such that r < f <C |6|. 
r+ is defined by j(r+) = J+. We have also approximated (1 — -j-)'^^ = e.~^ . The factor 
of 1/L above is the normalization of the wavefunction which we will not need in this 
particular computation. The "~" in the last two equations reflects an ambiguity by a 
constant prefactor that depends on the moduli at inflnity (which we set to one in section 
7 by rescaling the coordinates). 

The integrals in (8.13) can be solved (in terms of T-functions) by expanding in ri/|6| 
or r2/J+ yielding an answer in terms of a power-series in r. The lowest order term in 
the series is a constant so we find {hi{r)) ~ Oj + 6jr" with a > 0. This implies that the 
harmonics in the metric have the same small r behaviour so for r ~ 0, near the scaling 
point, we can evaluate the behaviour of the five dimensional metric in the decoupling limit 
(see e.g. [20, eqn. (2.15)]) and find the metric does not develop a throat. 

Cutting the throat 

Finally we would like to translate our insight above, that quantum mechanically one can- 
not reach the scaling point, and hence no infinitely deep throat develops, into a rough 
quantitative estimate of a mass gap in the dual CFT. 

We will do this in the particular example of the D6D6D0 three center solution in- 
troduced in section 6, which is scaling when N > 1/2. In section 6 we considered the 
non-scaling version of these solutions where N < 1/2 and we considered multiple DO cen- 
ters; here we work with a single DO center and scaling charges, but we will use the notation 
of that section. Note that, from eqns (6.1)-(6.2), this system is actually a particular in- 
stance of the general construction of section 7 since a = and a = b (in the notation of 
section 7). Thus our computation of the harmonic functions in the previous subsection 
applies to this system. 

The calculation will proceed in two steps. First, we want to estimate at which scale 

the quantum smearing cuts off the naive infinite throat. Second, we will translate this scale 

into a mass gap (in the dual CFT) by analysing a scalar field on a toy model geometry 

with a throat cut off at this scale. The scale at which we expect a deviation from the naive 

infinite throat is of the order of the minimum expected inter-center coordinate distance. In 

the case of the D6D6D0 scaling solution j+ = ^ and, furthermore, the angular momentum 

is related to the inter-center distances by (see the constraint eqns. (6.1) with /3 = | for 

asymptotic AdS space [20]) 

8jN N 

rgg = 8j , ro = jj-^ ^8jj, (8.17) 

where, in the second expression, we have made a large I approximation. As these expres- 
sions are linear in j we see that their expectation value in the n = state is directly given 
in terms of {j)oo = 1- So we find that quantum mechanically one expects the throat to be 
cut of at a scale of order e ~ -j > 1/2. 

A nice check on this estimate is to determine the charge dependence of the constant 
term in {h{r)) in the r ~ limit since, by a small r expansion, 

h{r) = —^ = - + 0{r) (8.18) 
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Note that this leading, r-independent, term feeds directly into the metric at small r so is 
a very relevant physical quantity to compute. A careful computation of this leading term, 
taking into account the normalization L in eqn. (8.13), yields 

(/^iW)~^ + 0(r) (8.19) 



with 7 a small number of order one. For the D6D6D0 system this gives (e^-*^) := (r^^ ) ~ -^ 

which confirms our previous computation. This shows that the estimate of e is relatively 

robust and does not depend very strongly upon which particular quantum expectation 

value we use to compute it. Moreover, as the metric and solutions are defined via the 

harmonics this expectation value very directly relates to the quantum expectation value of 

the metric. 

Now that we have understood the charge dependence of e we wish to translate this into 

a mass gap in the dual CFT. The computation of the mass gap in terms of a scalar wave 

equation on a capped-throat geometry is somewhat technical and has thus been relegated 

to appendix C. Here we will quote the final result for the mass gap A(Lo + -^o) in terms 

of e: ^ ^ 

A(Lo + Lo)~->-. (8.20) 

c Zc 

with c = p^ = 61, the central charge of the dual CFT. In the regime where A'^ ~ 1/2 (so the 

inequality above is saturated) this matches the expectation from the long string picture 

that the lowest energy excitation in the CFT is of order 1/c (see e.g. [19]). Whether the 

different scaling in the Cardy regime, N ^ I, reflects new physics of these solutions or is 

an artifact of our toy model geometry (see appendix C) would be interesting to explore. 

It would also be interesting to find an interpretation of the dependence of the mass 

gap on the parameter e ~ N/I r^ -^ which seems closely related to h = (Lq — c/24)/c. 

The latter played an interesting role as an order parameter in phase transitions in the dual 

CFT, as discussed in [20]. It was suggested there that this phase transition is due to a 

large number of winding modes being turned on. Although this is speculative, it might, 

itself, hint at a "long string" picture for the CFT at small h. It is interesting that the 

result (8.20) might also hint at such a picture. In any case these results and speculations 

only further highlight the importance of better understanding the dual M = (0, 4) CFT. 

9. Conclusions 

In this paper we have studied spaces of solutions of the field equations of A^ = 2 (A^ =1) 
supergravity in four (five) dimensions that describe multicentered black hole configurations. 
For suitable choices of the charges of the individual centers the resultant geometry is 
completely smooth so our setup also applies to spaces of smooth solutions. We have 
determined the appropriate symplectic form on these solution spaces and used that to 
quantize them and determine both the number of states and the explicit expressions for 
the wave- functions in several examples. 

The solutions we consider are all stationary, but in general not static, and, as such, 
they carry non-trivial angular momentum. Because of this momentum we did not need 
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to include small velocities for the centers in order to obtain a non-degenerate phase space 
(as in e.g. [34]). In other words, the space of stationary solution really is a phase space 
and not a configuration space. However, for some choices of the charges, we do find that 
the phase space is degenerate. This happens, for instance, when all centers carry DO and 
D4 brane charge so that all inner products {Tp,Tq) are identically zero (note that in this 
case there is also no angular momentum). In situations like these one could try to include 
small velocities for the centers in order to arrive at a well-defined phase space. It is clearly 
an interesting question whether this modified system will give rise to BPS states upon 
quantization. Superficially, the momenta increase the energy while leaving the charges 
invariant, and they therefore violate the BPS condition. However, if the Hilbert space 
has a continuous spectrum of momenta it is possible that there is a BPS bound state 
at zero momentum in the spectrum. This is difficult to analyze in general, but in our 
case we don't expect this to happen, at least not in asymptotically AdS spaces since AdS 
effectively provides a box and will therefore put an IR cutoff on the admissible momenta. 
It is crucial, for this argument, to be in the RR sector of the CFT. As was argued in 
[20], time independent solutions in the RR sector which include a purely fluxed D6 and D6 
center can be re-interpreted as possibly time-dependent solutions in the NS sector which no 
longer include the purely fluxed D6 and D6 centers. Therefore it is quite likely that there 
are interesting time-dependent BPS solutions in the NS-sector. We believe that the index 
computation of e.g. [59], where multicentered black hole configurations including velocities 
were considered, can also be reinterpreted along these lines and it would be interesting to 
work this out in more detail. Thus our proposal is that solution spaces with a degenerate 
symplectic form should not be thought of as describing proper BPS bound states. 

This immediately leads to another issue; namely, we know that, for example, N DO- 
branes can form a marginal bound state [60, 61], but the symplectic form for such a 
configuration (for example in the presence of a D4-brane) would vanish identically. This 
clearly conflicts with the statements of the preceding paragraph. 

We would like to argue that the resolution of this inconsistency lies in the fact that 
the marginal bound state of DO-branes cannot be understood purely from a low-velocity 
expansion in the Coulomb branch. Rather, the presence of the non-abelian degrees of 
freedom is essential for the bound state to exist. In other words, the marginal bound state 
has significant support in the Higgs branch. This is supported both by the analysis of [61] 
as well as by the size of the bound state (see e.g. [62]). Again, it would be interesting 
to explore this further. In this paper, we have taken the point of view that a solution 
containing a marginal bound state of e.g. N DO-branes should be counted separately 
from a similar solution where the marginal bound state has been replaced by N individual 
DO-branes, and we obtained satisfactory results. 

One of the main results of our paper is that we have explicitly shown how scaling 
solutions, which classically give rise to arbitrarily deep throats in space-time, are modified 
in the quantum theory. The deep throats are capped off in such a way that the spectrum of 
the dual CFT remains gapped, and we have estimated the size of the gap to be of order 1/c. 
This resonates well with the "fractionation" idea described in [19]. Here, the results suggest 
the possible existence of a long string picture, whose length is p^ times as long as that of 
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the original string. In any case, it is striking that our results show that geometries that 
are perfectly smooth and can have arbitrarily low curvature can nevertheless be ruled out 
as acceptable classical geometries. The reason, in our case, is that although the solutions 
develop a very deep throat in space-time, they are only responsible for a tiny fraction of 
the total phase space, somewhat analogously to what was observed in [18] . There is a nice 
qualitative picture one can use to understand this. The scaling solutions require localization 
of the centers within an arbitrarily small coordinate distance. If the centers did not warp 
the geometry, standard quantum mechanics (the uncertainty principle) would tell us this 
is not possible. The back-reaction of the centers, however, stretches this infinitesimal 
coordinate distance to macroscopic sizes. The novel physics emerges from the realization 
that the symplectic form does not seem to feel this "stretching" so the phase space volume 
stays small. 

To our knowledge, this is one of the first explicit examples where we can see that 
quantum gravitational effects are much stronger than what one would expect based on 
semiclassical reasoning. It is tempting to speculate that a similar phenomenon may take 
place at the horizon of a black hole and this is one of the key ingredients in resolving the 
black hole information paradox. 

The simplest scaling solutions are described by quiver quantum mechanical systems 
where the quiver has a closed loop. What we have studied is the quantization of the 
Coulomb branch of the theory. One could imagine that down the throat, where the branes 
approach each other very closely, a Higgs branch description is more appropriate, since 
the open string stretching between the branes become very light. However, in examples 
one finds that the number of states in the Higgs branch is much larger and can have 
exponential growth [15]. In these cases the Coulomb and Higgs branches are therefore 
either completely decoupled, or alternatively, only part of the Higgs branch can be seen on 
the Coulomb branch. This point clearly also deserves a better understanding. 

For a multicentered black hole configuration the quantization we are doing effectively 
counts the number of configurational degrees of freedom, but it clearly ignores any en- 
tropy the individual black holes might have. Therefore the resulting states do not have 
a straightforward interpretation as states in the dual CFT. The quantization of spaces of 
smooth solutions, on the other hand, does lead to BPS states which should have a dual 
description as BPS states in the CFT as well. At the same time, multicentered black hole 
configurations should also be describable by suitable ensembles of states in the CFT. One 
may therefore wonder to which multicentered black hole configuration a state obtained 
from quantizing a given space of smooth solution belongs. This is in particular of relevance 
for the "fuzzball" proposal [10] (see also review [19, 30, 63]) which, in its most optimistic 
form, would imply that all states that make up multicentered black holes can be obtained 
by quantizing suitable spaces of smooth supergravity solutions. 

In [15] it was proposed that the Hilbert space of states of a given total charge can be 
partitioned in sub-spaces each corresponding to a given attr actor flow tree. Furthermore 
we expect the states on a solution space of smooth centers to be elements of this Hilbert 
space. This seems to imply that such states should be related to a unique attractor tree. 
This is especially non-trivial in the case of smooth scaling solutions. Most naively one 
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could propose the following relationship between multicentered black hole attractor trees 
and smooth "multi-scaling" solutions. Let us assume that for a given set of charges there 
is only one way (in a given connected component of solution space) to cluster the charges 
such that the individual clusters form "local scaling solutions". By local scaling solutions 
we mean that in some regime in solution space all the charges in a given cluster can 
be brought arbitrarily close together to form an infinitely deep throat. Since there are 
multiple clusters there is also a regime in solution space where the charges in each cluster 
are brought together to form several separated throats. So our assumption is essentially 
that within a given connected component of solution space it is not possible to re-arrange 
the charges into a different set of clusters that still enjoy this property. Let us explore 
this possibility briefly as it may have interesting consequences. For instance, if this is true 
then it would be natural that the states arising from the quantization of such a connected 
component correspond to microstates of the multicentered black hole composed of centers 
each carrying a charge corresponding to one of the clusters that can form a scaling throat. 

The split attractor flow picture might then be useful in classifying this correspondence 
(between multi-black hole configurations and smooth multi-scaling solutions). Recall that 
all non-scaling multicenter solutions, smooth or not, are described by a suitable attractor 
flow tree of which the existence is necessary and (conjectured to be) sufficient for the 
existence of the multicenter solution. For the corresponding multi-scaling solutions the 
existence of the same attractor flow tree is necessary as well (although known not to be 
sufficient). In particular, when some centers locally form a scaling solution, which can 
have an arbitrarily deep throat, they cannot be distinguished individually in the attractor 
flow tree (because no amount of moduli tuning will force the centers apart and this is 
the basis of the attractor tree construction). The tree is thus conjectured [15] to end on 
the attractor point corresponding to the sum of the charges that locally form a scaling 
solution. The existence of the tree is necessary because when all the clusters have turned 
into infinitely deep throats the spacetime geometry away form an individual cluster is very 
well approximated by the geometry of a multi-black-hole solution with the same charges 
so if the latter has pathologies then so will the former. 

A tentative proposal is therefore that a given smooth solution belongs to the multi- 
centered black hole configuration which has the same attractor flow tree. The attractor 
tree either describes a multicentered black hole configuration with one black hole for each 
endpoint of the tree or it can be "resolved" to a similar configuration with each black hole 
replaced by a "local" smooth scaling solution. In particular, the smooth solutions that 
belong to a single black hole are those where all the centers can simultaneously disappear 
down a single, arbitrarily deep, throat; smooth solutions that belong to a two-centered 
black hole are those that can disappear down two throats; etc. So roughly each group 
of smooth centers that locally form a scaling solution belongs to a single black hole con- 
stituent. We should emphasize that while this is an appealing picture, very little is yet 
known about the structure of scaling solutions with more than three centers and this would 
be a prerequisite to making any such proposal precise. 

To test this idea and to see whether one can possibly account for the entropy of a 
single black hole in this way, we need to find scaling solutions with many centers, quantize 
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the corresponding space of solutions, and count the resulting number of states. The most 
promising example to study is the one consisting of D6D6 and many DO's studied in 
section 6. In section 6 we counted the number of solutions in the non-scaling regime 
and found perfect agreement with the wall-crossing formula. These solutions can also be 
arranged in such a way that they are scaling, but we have not yet completed a computation 
of the number of states one can obtain in this way (see also footnote 16). We hope to report 
on this in the near future. 

It is also interesting to note that our solution space quantization is somewhat remi- 
niscent of the "near horizon probe" proposal of [56, 64, 58] and hence might lead to some 
insights in the exploration of the link between this proposal and the fuzzball picture [65, 66] . 

Finally, there are several other loose ends and open problems to explore. These include 
a proper understanding of the dual CFT and the right way to describe the CFT duals of the 
multicentered black hole configurations. We hope that the results in this paper, such as the 
computation of the gap, will allow us to make some progress in this direction. We would 
also like to have a better understanding of the role of the fermions in the quantization, 
and like to know whether the solution space always comes out to be Kahler or even toric 
(which is supported by the agreement with and the derivation of the wall-crossing formula 
in [15]). It would also, perhaps, be interesting to compute the symplectic form directly in 
supergravity and understand why terms that do not come from the gauge field cancel (see 
discussion in the beginning of section 4). We leave these, and several other issues, to a 
future study. 



Acknowledgements 

The authors would like to thank Frederik Denef for initial collaboration and discussion 
during later stages of the project. We also gratefully acknowledge useful discussions with 
I. Bena, M. Berkooz, T. Levi, S. Mathur, G. Moore, K. Papadodimas, B. Pioline, and 
S. Raju. 

DVdB is an Aspirant of the FWO Vlaanderen and is partially supported by the Euro- 
pean Community's Human Potential Programme under contract MRTN-CT-2004-005104 
'Constituents, fundamental forces and symmetries of the universe', by the FWO - Vlaan- 
deren, project G. 0235. 05 and by the Federal Office for Scientific, Technical and Cultural 
Affairs through the 'Interuniversity Attraction Poles Programme - Belgian Science Policy' 
P6/11-P. The work of JdB, SES and IM is supported financially by the Foundation of 
Fundamental Research on Matter (FOM). 



A. The Three Center Solution Space 



In this appendix we will analyze some properties of the moduli space of three-center so- 
lutions. Our starting point will be the set of equations (2.6) which we will rewrite as 
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y 


z 


c 


a 



follows 

a b 

- - - = C2-C1 

C3 -C2 

C1-C3. (A.l) 

Here a,b,c represent the inner products {Ta,Tf,), x,y,z are the lengths of the three sides 
of the triangle spanned by x^, and C2 — ci = {h,Ti) etc. The constants Ca are not uniquely 
fixed, as we shift them by a fixed amount without modifying the above equations. Still, 
expressing things in terms of Ca allows for a somewhat more symmetric treatment. 

The first important remark is that up to an 50(3) rotation, x, y, z uniquely determine 
the solution. In other words, the quotient of the solution space by 50(3) is precisely the 
set of solutions x,y, z of (A.l). 

Second, we should keep in mind that x, y, z are the sides of a triangle, i.e. they 
should be nonnegative numbers that satisfy the triangle inequality x + y > z and its cyclic 
permutations. 

In our discussion of the solution space quantization, the size of angular momentum will 
play an important role, as we will use it as a coordinate on the solution space. In terms of 
the variables used in (A.l), the angular momentum is given by (see (2.11)) 

J^ = -- (X^(C2 - Cl)(ci - C3) + y^(c3 - C2)(C2 - Cl) + Z^(ci - C3)(c3 - C2)) . (A.2) 

We would in particular like to know whether | J| is a good single- valued coordinate on the 
solution space and what range of values it takes. 

It is easy to write down the general solution to (A.l) in terms of a single free parameter 
A: 

?/=V-^' ^ = V^- (A-3) 



A — Cl A — C2 A — C3 

This is the general solution if a, b, c are not equal to zero. If all three are zero, or two out 
of three are zero, there are either no solutions to (A.l) or the space of solutions is at least 
two-dimensional. In either case the symplectic form becomes degenerate and most likely 
these solution spaces do not give rise to BPS states^^. Finally if one of a, b, c is zero, say 
a = 0, then either there are no solutions or one finds a fixed value for y, z from (A.l), while 
X is not constrained by (A.l). However, x is constrained by the triangle inequalities so 
that the solution space becomes 

a = 0, & / 0, c / 0, =^ y, z fixed, |y — ^| < x < y + z. (A.4) 

We now continue with the case where a, b, c are not equal to zero so that the solutions 
are of the form (A. 3). We again need to distinguish a few cases. The most degenerate 



^^In order to obtain a non-degenerate symplectic form we would have to add additional degrees of freedom 
to the theory, which would act like momenta for some of the coordinates. A nonzero momentum will typically 
increase the energy while the charges remain fixed, violating the BPS condition. Therefore we expect no 
BPS states in this case. 
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case is when ci = C2 = C3. Then either the moduh space is empty or one-dimensional, but 
in the latter case the angular momentum vanishes identically everywhere on the solution 
space and thus the symplectic form is trivially degenerate. 

The next case is a, b, c nonzero and two of the Cj equal to each other. Using the 
permutation symmetry of (A.l) and the possibility to simultaneously change the signs of 
a,b,c,Ci, X, we can distinguish three different cases: (i) C3 > ci = C2 and a,b,c > 0, (ii) 
ci = C2 > C3 and a,b,c> 0, and (iii) C3 > ci = C2, a, 6 > and c < 0. Positivity of x, y, z 
requires that A G /i = (03,00) for cases (i),(ii) and A S /i = (01,03) in case (iii). Next we 
denote by I2 the set of solutions of the triangle inequalities 

a + b c \a-b\ , . ^, 



A — ci A — C3 A — ci 

It is easy (though somewhat tedious) to see that h n I2 is either empty, an interval of the 
form [A- , A+] , an interval of the form [A_ , 00) , an interval of the form (ci , A+] , or an interval 
(ci, 00). The endpoints A+ and A_ always correspond to a point where a triangle inequality 
is saturated. The interval extends all the way to infinity only if"^*^ a + b > c > \a — b\, i.e. 
when a, b, c satisfy triangle inequalities, which can only happen in case (i) and (ii). In these 
cases there is a scaling throat with x,y,z -^ 0. The interval starts at ci only if we are in 
case (ii) or (iii) and a = b, and in this case the solution space includes configurations where 
a center can move off to infinity. 

From the point of view of angular momentum, the case where one of the centers moves 
away to infinity (e.g. x, y — > 00) can be viewed as a case where the triangle inequalities 
X + z > y and y + z > x are both saturated. Therefore, in all cases we have analyzed 
so far, the solution space contained just a single connected component described by single 
interval of possible values of A, and at the endpoints of the interval either one has a 
scaling solution with vanishing angular momentum, or a solution that saturates at least 
one triangle inequality. Whenever this happens, we always find that 

|J|2 = i(±a + ±6+±c)2 (A.6) 

for suitable choices of the signs, as can be seen easily e.g. from (2.7). 

It remains the analyze the generic case with all Cj different from each other. Up to an 
overall sign flip and a permutation, there are two cases, which are (iv) ci < C2 < C3 and 
a, 6, c > and (v) ci < C2 < C3 and a, 6 > 0, c < 0. Positivity of x, y, z in case (iv) implies 
A > C3 and implies C2 < A < C3 in case (v) . The main problem is to analyze the triangle 
inequalities. They can be analyzed qualitatively by sketching x + y — z, x — y + z and 
—X + 2/ + 2; as a function of A. We know that each of these functions can have most two 
zeroes as a function of A, and we know its behavior near the three poles at A = ci, 02,03. 
We will skip the details, but one finds that the moduli space consists of at most two 
components, each of which corresponds to a certain interval of possible values of A. At the 
boundaries of each interval a triangle inequahty is saturated. Notice that in case (iv) one 



* Actually, this can possibly also happen when a-\-h — core— |a— &| 
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of the components can be of the form [A_, oo). This is possible whenever a, b, c themselves 
satisfy triangle inequalities. If this happens, at A = oo there is a scaling solution. 

To summarize, the solution space in all cases consists of at most two components, 
corresponding to two intervals of possible values of A. At the endpoints of the interval 
some triangle inequality is saturated. This can include configurations where one of the 
centers moves off to infinity (cases (ii) and (iii) above, with a = b), and scaling solutions 
where A — > oo (cases (i), (ii) and (iv) with a,b,c obeying triangle inequalities). 

Finally, we would like to show that J^ is a good coordinate on each component of the 
moduli space of solutions to (A.l). In order to do so, we compute dJ^/dX. According to 
(A. 3), dx/dX = —'J? I a and similarly for y,z. If we differentiate (A. 2), use these relations, 
and finally replace q by the left hand side of the original equations (A.l), we obtain 

^d/^ ^ i!i! /« _ ^^ /c _ «\ _^ r^ 1^^ _ c\ /a _ 6\ z^ /c _ a\ /6 _ c\ 
dX a \x yj^z x) b \y z J \x yj c \z xJ \y zj 

We rewrite this as 

dJ^ 2 f ^i \ 4n2no — n? 

—zabcxyz—— = n^a + nia + n2 = uq \ a-\ H (A. 8) 

dX \ 2nQj AriQ 

with no,ni,n2 certain a-independent polynomials. The right hand side of (A. 8) is positive 
if no and 4n2no — n\ are positive. By explicit computation we find 



„^ _,.-., ^y,.2 .,2 .2^^^, cV 



( 2i , ^y ( 2 2 2\\ I 

[z b+Y^^ -y -^ )) + 



4z^ 
4n2no -nj = b'^c^x^{bz - cyfO (A.9) 

where 

6 = {x + y + z){x + y - z){x - y + z){-x + y + z). (A.IO) 

Since > if all triangle inequalities are satisfied, we have indeed shown that J^ is a 
monotonous function of A and that J^ is a good coordinate on each component of the 
solution space. 

B. Aspects of Toric Geometry 

In this appendix we review some techniques in toric geometry that we use throughout the 
paper. By construction we start in our description of solution space in the main text from 
a symplectic point of view. It is however more convenient for geometrical quantization 
to have a Kahler description, which can always be made in the case of a symplectic toric 
manifold. The main formulas in this appendix are thus the expressions (B.5,B.6) for the 
complex coordinates and Kahler potential in terms of the symplectic coordinates on a 
symplectic toric manifold. Before giving these formulas we review some of the basics of 
symplectic toric manifolds and symplectic toric orbifolds. 
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Polytopes 

As is customary we will refer to the convex hull of a finite number of points in M" as a 
polytope. The boundary of such a polytope is itself the union of various lower dimensional 
polytopes that are called faces. In particular a zero-dimensional face is called a vertex, a 
one-dimensional face an edge and a n— 1-dimensional face a facet. Note that we can view any 
polytope as the intersection of a number of affine half spaces in M". A polytope P can thus 
be uniquely characterized by a set of inequalities, namely x G P iff Vi = 1, . . . , ^^(facets) 

(q,x) > Aj <^ X^CjjXj > Aj. (B.l) 

j 

Given a polytope we will call the set Ci £ Z"', given by the inward pointing normals to the 
various facets, the normal fan. 

An n-dimensional polytope is called a Delzant polytope if it satisfies the following three 
conditions 

• simplicity: in each vertex exactly n edges meet, 

• rationality: each of the n edges that meet at the vertex p is of the form p + tui with 
t G M+ and m G Z", 

• smoothness: for each vertex the Ui form a Z-basis of Z"". 

The polytope is called rational instead of Delzant if we replace in the third condition the 
requirement of a Z-basis by that of a Q-basis. 

Symplectic Toric Manifolds 

Before giving the precise technical definition of a symplectic toric manifold, let us first 
sketch the idea. Roughly speaking a toric manifold is a T" fibration over a given n- 
dimensional polytope, such that at each facet a single U{1) inside the T" shrinks to zero 
size. On the intersections of the different facets multiple [/(l)'s collapse, e.g. at the vertices 
all circles have shrunk. On the interior of the polytope the toric manifold is simply of the 
form P° X T" and the full toric manifold is the compactification of this space. On the 
interior there is thus a standard set of coordinates of the form (xi,9i) with Xj G P° and 
0j G T and the manifold comes with a standard symplectic form 0, = ^^ dxi A d9i. It is of 
course rather non-trivial that this manifold can be smoothly compactified, but when the 
polytope is Delzant this is the case. Let us now state the above ideas more precisely. 

A symplectic toric manifold is a compact connected 2 n-dimensional symplectic man- 
ifold (M,0,) that allows an effective Hamiltonian action of an n-dimensional torus T". 
Remember that the action of a Lie group on a symplectic manifold is called Hamiltonian 
if there exists a moment map fi from the manifold to the dual Lie algebra that satisfies 

d{fi{p),X)=n{;X), (B.2) 

with p G M, X a generator of the Lie algebra and X the corresponding vectorfield. Fur- 
thermore the moment map should be equivariant with respect to the group action, i.e. 
f^idip)) = Ad* o fi{p), with Ad* the coadjoint representation. 
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By a theorem of Delzant [67] every symplectic toric manifold is uniquely characterized 
by a Delzant polytope. Given a symplectic toric manifold the corresponding polytope is 
given by the image of the moment map. To conversely reconstruct the manifold from 
the polytope is slightly more involved and relies on the technique of symplectic reduction, 
we refer readers interested in further details to e.g. [68]. Note that the normal fan to the 
polytope can be interpreted as a fan, which is used to characterize toric varieties in algebraic 
geometry, see e.g. [69] for a nice introduction. This can be useful to identify a symplectic 
manifold given by a polytope and furthermore provides an embedding in projective space. 

Kahler toric manifolds 

What will be of use in this paper is that Delzant 's construction also associates a set of 
canonical complex coordinates to every symplectic toric manifold, effectively implying that 
every closed symplectic toric manifold is actually a Kahler manifold. As throughout the 
paper we will make use of the explicit construction of these complex coordinates in terms of 
the symplectic coordinates (xj, 9i), we will detail the general procedure here, be it without 
proofs or motivation. Those can be found in references [45, 46]. 

As mentioned above (B.l) any polytope P is characterized by a set of inequalities. 
Given this combinatorial data of the polytope one can define associated functions 

which are everywhere positive on P. Using these functions one can define a 'potential' as 
follows 

9{x) = lY.^i{x)logk{x). (B.4) 

i 

In case the polytope is Delzant, it is shown in [45] that this potential defines good complex 
coordinates on the toric manifold as follows 

Zi = exp ( -^g{x) +i9i\ . (B.5) 

Furthermore a Kahler potential for the corresponding Kahler metric $!(•, J-) is given by 

/C = J^Aaog/.(x) + ?oo. (B.6) 



It follows from the construction [45, 46] that /C is the Legendre transform of g, i.e. fC{z) 

dx' 



—,x — g{x). This can be used to derive that 



de.a.<>;K = exp(i:g)det^, (B.7) 



which will be a useful formula in the bulk of the paper. 
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Toric orbifolds 

As we also consider quotients of symplectic toric manifolds by a permutation group in 
this paper, it will be necessary to introduce the generalization of the above construction 
of complex coordinates to that of symplectic toric orbifolds. As in the manifold case, a 
symplectic toric orhifold is a 2n-dimensional symplectic orbifold that allows a Hamiltonian 
T" action. As was shown in [54] such symplectic toric orbifolds are in one to one corre- 
spondence to labeled rational polytopes. Such a labeled rational polytope is nothing but 
a rational polytope with a natural number attached to each facet. The label mi denotes 
that the i'th facet is a Z^^ singularity. Again the explicit construction of the toric orbifold 
from the labeled polytope is rather involved and we refer those who are interested to [54]. 
The labeled polytope corresponding to the quotient of a symplectic toric manifold by a 
group respecting the torus action, is however easy to find. It is given by the quotient of 
the original polytope and attaching a label m to each facet that is a Z^ fixed point under 
the group action. 

Given a labeled rational polytope one can construct complex coordinates on the toric 
orbifold in a way similar to the manifold case. The functions U from (B.3) are generalized 
to [54, 55] 

li{x) = "T-j ^ CijXj - Aj , ^oo = ^ ruiCijXj , (B.8) 

\ J J i,j 

where rrii is the label attached to the facet orthogonal to the vector Cj. The construction 
of the complex coordinates and the kahler potential from these functions then carries on 
analogously to (B.5,B.6). 

C. Gravitational Throats and CFT Mass Gaps 

Having an infinitely deep throat in AdS spaces seems paradoxal as it suggests a contin- 
uous spectrum on the CFT side, as one can make arbitrarily small energy excitations by 
localizing them deep enough in the throat. In this appendix we would like to calculate 
the correspondence between the size of the mass gap and a smoothly capped off throat 
in the bulk. We are going to approximate the finite throat by a toy model metric with a 
throat that is crudely cut off at a scale e and try to solve the scalar wave equation in this 
background. Due to the technical difficulty of the full problem we resort to matching the 
far (r ^ e) and near (r <C e) region solutions. By doing so we are able to relate the inverse 
of the depth of the throat to the mass gap on the CFT. The result we find agrees with 
what one expects for CFT's with a long string picture. We use this result in section 8 to 
relate the cutoff scale e we calculate there to the size of a mass gap. 

To get a reasonable guess to what the capped off geometry would be, we use that far 
away from the tip of the throat the geometry looks like the one of a D4D2D0 BTZ black 
hole. So our starting point is the following metric (see (4.1) in [20]) 



r 



Ir + C , ,n ^^ndr^ 



ds^ = --dtdil: + -—-—dtl}'^ + U^^^ (C.l) 

U A U r2 
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with C = 2„3 some constant determined by the charges. Let us assume this metric to be 
a good approximation for r > e, whereas we take 

2tt2 tt2 

ds^ = --dtdip + -^^dip'^ + -^dr^ (C.2) 

U e^ e^ 

vahd for r < e. We find the fohowing radial equations for a free scalar field in these two 
'sub-geometries', in the outer region 

-a.^5. + m^ - ^, —j <P = (C.3) 

and in the inner region 

-drjj^dr+m^ -^ _J0 = O. (C.4) 

Let us solve these equations. The first one has two solutions but only one is normalizable 
at infinity. Define A via 

1 _ A2 = -m'^U^ (C.5) 

4 



and also 



C^-i|V^('^ + 4fe) (C.6) 



and the variable 

z = 2uiVUC/r (C.7) 

then the field equation for (p becomes 

(^' + (-^ + ; + i^))'^ = (C.8) 

and the solution in the outer region in terms of Whittaker functions is (see e.g. eq (84) in 
[70]) 

(/.(z) = A%a(^) + M„5,a(-^). (C.9) 

These two terms are proportional to each other but this answer is the linear combination 
which yields a real answer. The large z behavior is 

Next we turn to the cap region. Here, the field equation is usually solved in terms 
of parabolic cylinder functions. We need the right linear combination which vanishes as 
r — > to be smooth there. Define 



2 / '^ 

em , / e'^ 



2fe, r = x\-—: (C.ll) 



2U ' V 2u;C/3 
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,, . , -^^ + c.c. (C.12) 



then the field equation is {dy + ^ — o)(/) = and then we find the following large x behavior 
of the solution (notice that r = e is at large x) 

i n3/8n-l/4-ia/2 

^^exp(-..V4)x-^/^ ^-^^ ^,,3 ^,. 

^4 "T 2 ^ 

Notice that the coefficients are important, only for this particular coefficient (and including 
the c.c.) this function has the asymptotics of the regular solution. 

In order to match the solutions (C.9) and (C.12) at r = e, both the fields and their 
first derivatives should match. A priori this need not be possible since we already have 
unique solutions up to overall normalization. Since the overall normalization is not fixed, 
the condition that we need to impose is 

dz (J)'{z) , _ ^dx (p'jx) ^ ^ ^^ 

dr (j){z) ^~'^ dr (f){x) ^~'^' 
In the small e-limit, the derivatives of the exponentials give the largest contributions to (p' , 
thus we only need to differentiate those. 

To proceed, we define e^f to be the phase of l/r(l/2 + A — ^), and e*'^ to be the phase 
of 

(_l)3/82-l/4-W2gaV4^ 

r(| + f) • ^ • ^ 

Then we find 

z 1 X 

(pouter ^ COs{— - {^/i)log{z/i)+ p), (/)inncr ~ ^ C0S( + a log(x) + Cj) (C.15) 

Zl \ X 4: 

up to irrelevant overall normalizations. We now evaluate the matching condition keeping 
only the leading terms for small e, i.e. we only differentiate z/2i and — a;^/4. After doing 
this, various factors of e and ui happily cancel and we are left with the matching condition 
(evaluated at r = e) 

C/5/2 x^ Z 

-^^ tan(-— + alog(x) + a) = tan(- - (^/i) log(zA) + p). (C.16) 

As we increase loj, but keeping oo small enough so the matching strategy remains sensible, 
the ffist tangent seems to oscillate most rapidly (as long as U^ > 2\/UC, otherwise the 
other tangent seems to win). Because of this oscillatory nature we get a gapped spectrum. 
A crude estimate for the gap can therefore be made by looking at the values of a; for which 
the ffi'st tangent makes a n period. Since x^/4, evaluated at r = e, is ujU'^/{2e), we find 
the following rough estimate: 

A^~^. (C.17) 

Since the eigenvalue of dt is like that of Lq + Lq, without any further factors, we finally 
conclude that 

A(Lo + Lo)~|j. (C.18) 

In case e is of order 1 the gap in the conformal weights scales like 1/U^ which is precisely 
1/c, with c the central charge of the dual CFT, as one would get from long string fraction- 
ation. This suggests that long strings might play an important role in the physics of four 
dimensional black holes and the dual A^ = (0, 4) SCFT. 
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